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ABSTRACT AND SUMMARY 


This report contains a number of preliminary studies on the 
application of a Rotating Gravity Gradiometer (RGG) system on 
board a Lunar Polar Orbiter (LPO) for the measurement of the 
Lunar gravity field. An RGG system would be desirable for the 
LPO mission since it would significantly reduce the total mission 
weight and cost budgeted for the gravity portion of the science 
by eliminating the need for a Relay Satellite and by decreasing 
the tracking time needed for full Lunar coverage from months of 
continuous tracking to daily data dumps. However, there were a 
number of questions on the compatibility of the present airborne 
RGG design with the LPO mission. The various studies in this 
report address those questions. A data collection simulation 
study shows that a gradiometer will give significantly better 
gravity data than a doppler tracking system for the altitudes 
under consideration for the LPO, that the present demonstrated 
sensitivity of the RGG is adequate for measurement of the Lunar 
gravity gradient field, and that a single RGG instrument will 
provide almost as much data for geophysical interpretation as an 
orthogonal three axis RGG system. An engineering study of the 
RGG sensor/LPO spacecraft interface characteristics shows that 
the RGG systems under consideration are compatible with the pre- 
sent models of the LPO spacecraft, can be placed anywhere on the 
spacecraft, do not require any special alignment or thermal envi- 
ronment, do not generate RF, magnetic or other disturbing fields 
and put no new constraints on spacecraft attitude, rate and con- 
trol system accuracies. The interface study also contains a 
physical and functional description of a typical RGG system as 
designed for a spacecraft. An error simulation study shows that 
the errors in spacecraft orbital position and spacecraft attitude 
will not significantly degrade the gravity gradient data during 
the data reduction process. 
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CONCLUSIONS 


The general conclusion of this multi-study effort is that a 
Rotating Gravity Gradiometer is not only compatible with a Lunar 
Polar Orbiter mission, but will provide science data of higher 
resolution and higher sensitivity than the planned Doppler track- 
ing system. It will also significantly reduce the total mission 
weight required for the gravity experiment by eliminating the 
need for a Relay Satellite and will decrease the DSN tracking 
time required for full Lunar coverage from nearly continuous 
tracking for months to daily data dumps. In Appendix A, a data 
collection simulation study shows that a gradiometer will give 
significantly better gravity data than a doppler tracking system 
for the altitudes under consideration for the LPO. The study 
also simulated the effects of sensor noise on the mapping of 
gravity gradient data and the integration of gradient data into 
gravity data and found that the present demonstrated sensitivity 
of the RGG is adequate for measurement of the Lunar gravity grad- 
ient field. Other portions of the data simulation study found 
that a single RGG instrument oriented along the orbital axis to 
measure the gravity gradients in the orbital plane will provide 
almost as much data for geophysical interpretation as an orthogo- 
nal three axis RGG system. The cross gradient output of the RGG 
should be integrated with the spacecraft velocity to obtain a 
contour map of the vertical gravity field, and then the higher 
resolution (but slightly distorted) principal gradient difference 
output of the RGG should be used to enhance the resolution of the 
data set for geophysical interpretation purposes. In Appendix B, 
an engineering study of the RGG sensor/LPO spacecraft interface 
characteristics shows that the RGG systems under consideration 
are compatible with the present models of the LPO spacecraft, can 
be placed anywhere on the spacecraft, do not require any special 
alignment or thermal environment, do not generate RF, magnetic or 
other disturbing fields, put no new constraints on spacecraft 
attitude and rates, and should not affect the spacecraft attitude 
control system if the RGG angular momentum is included in the 
spacecraft control equations. Appendix B of the report also con- 
tains a physical and functional description of a typical RGG sys- 
tem as designed for a spacecraft. Another simulation study (See 
Appendix C and D and Task 5) that includes the orbital motion and 
attitude variations of the spacecraft shows that the variation in 
spacecraft orbital altitude and position, and spacecraft attitude 
errors will not significantly degrade the gravity gradient data. 
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RECOMMENDATIONS 


Because of the compatability of the RGG system with the LPO 
spacecraft and mission, and because of the significant advantages 
to be gained by replacing the Relay Satellite and Doppler veloci- 
ty tracking system with an on board gravity gradiometer, it is 
recommended that serious consideration be given to the funding of 
an instrument development effort to produce a breadboard version 
of an RGG sensor and spacecraft compatible electronics, combined 
with sufficient testing and demonstration of the sensor to reduce 
the technical risk associated with a new instrument concept and 
to increase the confidence level in the LPO Science Board on the 
use gradiometer instruments as a means of measuring Lunar gravi- 
ty. 
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STUDY TASK RESULTS 


TASK 1 - MEASUREMENT SENSITIVITY AND REPEATABILITY 


The first task on the study was to evaluate the 
sensitivity and repeatability of gravity gradient 
using the existing Hughes Rotating Gravity Gradiometer 
trolled laboratory conditions. 


measurement 
measurement 
under con- 


Since the Hughes RGG was undergoing continuous testing, 
rework and improvement during the study period, this was a con- 
tinuing effort throughout the contract. 


Measurement Sensitivity 


A single Rotating Gravity Graoiometer measures the gravity 
gradient components in the plane of its rotation. In operation, 
the phase of the signal demodulation process is adjusted so that 
the zero phase of demodulation occurs when the rotor index is 
aligned with one of the sensor external mounting reference direc- 
tions. The demodulated RGG signal then contains two independent 
outputs, the sine (Sxy) or inphase output which contains a mea- 
sure of the difference in the two principal gravity gradient com- 
ponents in the sensor plane of rotation (Gxx-Gyy) , and the cosine 
(Cxy) or quadrature phase output which contains a measure of 
2 Gxy or twice the cross gradient component in the plane of rota- 
tion. 

If three orthogonal RGG sensors are used, then the five 
independent components of the gravity gradient tensor can be 
obtained from the six sensor outputs by combining the sensor out- 
puts as follows: 

Gxx = (Sxy - Szx)/3 

Gyy * (Syz - Sxy)/3 

Gzz = -Gxx -Gyy (not independent) 

Gxy = Cxy/2 

Gyz « Cyz/2 

Gzx * Czx/2 
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Since the Laplacian of the gravitational gradient tensor is zero, 

Gxx + Gyy + Gzz = 0 

and this is just the trace of the gravity gradient tensor, we 
find that the principal components of the tensor are not indepen- 
dent. 


When an analysis of the propagation of the sensor random 
noise error through this computational process is carried out, we 
find that the error in calculating the gravitational gradient 
tensor from the sensor outputs is slightly different for the two 
types of gradient components since we have six sensor outputs and 
only five independent gradient components. For the cross gradi- 
ent components, the measurement sensitivity is just half of the 
sensor random noise in the cosine output. For the principal com- 
ponents, the measurement sensitivity is 1.414/3 of the sensor 
random noise (since tnere are two independent sensor measurements 
combined to obtain eich principal gradient component) . This 
ratio is 1/2.12 or a little less than half of the sensor random 
noise in the sine output. 

Periodically during the contract, at times when the sensor 
was not being tested or reworked on our ongoing AFGL/ARPA Rotat- 
ing Gravity Gradiometer Development contract, we would carry out 
sensitivity and repeatability measurements with the sensor ori- 
ented so that its spin axis was in the vertical orientation. In 
this orientation, the effects of earth gravity on the sensor were 
a minimum. They are still not negligible however. At its pre- 
sent stage of development, the RGG still contains a large sensi- 
tivity to earth gravity due to its anisoelastic response. To 
correct this error sensitivity requires cutting metal in the sen- 
sor arm plates and this final trim is not planned for a number of 
months until after the final, high precision bearings are deliv- 
ered and installed and the sensor error coefficients remeasured 
on the new bearings. Despite the fact that there is still a 
large earth gravity sensitivity in the present configuration, we 
were able to obtain error sensitivity data that was almost at the 
final sensitivity goal of +1 Eotvos for a 10 sec integration 
time. (1 Eotvos is defined as~*l nanogal/cm where a gal (Galileo) 
is 1 cm/sec/sec. An equivalent unit for lunar mascon work would 
be 0.1 mgal/km or 0.1 milligal accuracy of measurement of lunar 
gravity due to a kilometer change in distance.) 

The plots shown in Figure 11) are the two outputs of the RGG 
in the laboratory with the spin axis vertical showing the random 
noise characteristics over a fifteen minute period. The sensi- 
tivity levels quoted are the sensitivity in calculating the grav- 
ity gradient component from the sensor output, which is roughly 
half of the unprocessed sensor output. 
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Repeatability 


1 


Repeatability of the RGG sensor is best demonstrated in the 
laboratory by moving a test mass close to the sensor, thus 
inserting a gravity gradient signal into the sensor to increase 
the sensor output, then later removing the mass, measuring the 
decrease and comparing its magnitude with the previously caused 
increase. It is not possible with the present laboratory setup 
to replace the mass accurately enough to repeat the test again 
with exactly the same gravity input, since a change in placement 
of as little as 100 microns at the 20 cm distances that separate 
the test mass from the sensor can cause a significant error in 
the gravity input. 

An example of tne repeatability of the RGG in tne laboratory 
is given in Figure {2}. The magnitude of the sensor output was 
first determined by taking a 100 sec average, then the test mass 
was placed in proximity to the sensor. After the signal rise had 
settled, the magnitude of the sensor output was measured again 
with another 100 sec average. The test mass was then removed, 
and after the fall time had settled, a tnird 100 sec average was 
taken. The signal rise was +417.3 Eotvos, while the signal fall 
2.5 minutes later was -416.7 Eotvos, or a repeatability to within 
0.6 Eotvos (within tne random noise). 



Figure 12} - Repeatability of RGG in Laboratory 
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TASK 2 - IDENTIFICATION AND TABULATION OF ERROR SOURCES 

Our work on Task 2 involved tne identification and tabula- 
tion of the error sources associated with the gravity gradient 
measurements of Task (1) . Those errors associated with the oper- 
ation of tne gradiometer in the Earth's gravitational field were 
then separately identified. By separating out the errors caused 
by operation in the Earth's gravity field, we can then make a 
better estimate of tne expected performance of the RGG sensor in 
the nearly free-fall LPO spacecraft environment. 

The details of the various sensor error sources are docu- 
mented in the ongoing contract reports [1,2] on the development 
of the Rotating Gravity Gradiometer. Most, it not all, of the 
error sources of concern are due to the effects, either direct or 
indirect, of the one g environment of the earth acting on the 
sensor. Of tne 23 separately identified error sources, 14 are 
directly or indirectly dependent upon the earth's gravity field. 
The only error sources that are independent of the earth's gravi- 
ty field are: thermal noise (which is the largest error source at 
1/3 the error budget) , speed control servo noise, data digitali- 
zation noise, scale factor or gain shift errors, sensor resonant 
frequency shift, and magi "stic field sensitivity. 

After some consideration and identification of the effect on 
the error coefficients of the earth's gravity field. Table 11) 
was prepared to show the expected contributions to the sensor 
errors of the various sensor error coefficients assuming the more 
benign environment of tne LPO spacecraft. 



TABLE 1 


OF RGG ERROR SOURCES 


Worst Case 
Airborne Prototype 
(Eotvos (§ 10 s) 

Estimated 
Lunar Orbiter 
(Eotvos @ 1 0 s) 

0. 30 

0. 30 

0. 26 

0. 26 

0.23 

0. 05 

0.21 

0. 21 

0.20 

0. 05 

0. 20 

0. 00 

0. 10 

0. 10 

0. 15 

0. 15 

0. 14 

0. 00 

0. 12 

0. 00 

0. 10 

0. 05 

0.64 

0.49 


T1732 


12 





TASK 3 - SCIENCE RETURN TRADE-OFF STUDY 

One of the first studies carried out in this program was a 
trc-.e-off study to determine the science return (both sensitivity 
and gravity harmonics (resolution]) as a function of orbital alt- 
itude. This was done by carrying out a computer simulation of 
what we could expect from a typical RGG system operating in on 
LPO spacecraft at various altitudes compatible with the LPO mis- 
sion. The most important concern was the type (format, direc- 
tionality and units) and quality (sensitivity and resolution) of 
the science return that we could expect from the total system. 
Was it adequate for scientific purposes? Was the data in a form 
that was readily utilized by other scientists not in the special- 
ization of gravity gradiometry? we also broadened the trade-off 
study to include the question: How did the data ^otaincd from a 
gravity gradiometer system compare in type and quality to that 
obtained from other methods of obtaining the same information 
(specifically that obtained from a Doppler velocity tracking sys- 
tem that was the original choice of the LPO mission analysts 
because of its previous utilization in mapping Lunar mascons)? 

The results of the study are presented in Appendix A - Grav- 
ity Gradient Mapping from the LPO - Science Return Simulation 
Study by Dr. Robert L. Forward, the Principal Investigator (with 
substantial assistance from Richard W. Lowe) . The paper has been 
accepted (without change) for publication in The Moon . The 
results of the simulation were very positive and surprising to 
everyone (including the highly biased PI) . • 

The simulations showed that gravity gradient data obtained 
using existing designs of airborne gravity gradient sensors would 
be of high resolution and high sensitivity and thus of the high 
quality needed for science interpretation. The simulations 
showed that certain components of the gravity gradient data could 
be plotted directly into gradient contour maps and used by geo- 
physicists for interpretation of subsurface structures, or alter- 
natively, other components of the gravity gradient could be easi- 
ly converted into more familiar gravity contour maps. Also, 
gravity contour maps obtained from gravity gradient data were 
found to have better resolution and signal-to-noiee than gravity 
contour maps obtained from Doppler tracking data. 

In the simulation studies, a frontside map of the Lunar mas- 
cons taken from a paper by Sjogren, Wimberly and wollenhaupt (3] 
were used to generate a point mass model. The point mass model 
was taken as the starting basis for the simulation since it could 
reproduce a gravity field that had a very close resemblance to 
the actual Lunar gravity field and yet had a number of high fre- 
quency features that could test the resolution and sensitivity 
requirements of any gravity measurement system. The results of 
the work are summarized by the abstract of the paper: 
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AESTFACT 


Simulations of the gravity data to be expected from 
a Lunar Polar Orbiter spacecraft utilizing either a Dop- 
pler velocity tracking system or a gravity gradiometer 
instrument system arc generated using a point mass model 
that gives an excellent representation of the types of 
gravity anomalies to be found on the moo . If the state 
of the art in instrumentation of both systems remain at 
the level of 1 mm/sec at 10 sec integration time for 
the Doppler velocity system accuracy and at + 1 Fotvos 
(1 Lotvos = 10 nrad/sec/sei, - 1 nqal/cm = 0.1 mgal/Km) at 
10 sec integration time for the gravity gradiometer sys- 
tem accuracy, inspection of the simulation plots show 
that a gravity gradiometer system will give science data 
with resolution that is nearly twice that of the Doppler 
system along with signal-to-noise ratio that is twice 
that of the Doppler velocity system at altitudes be' i 
100 Km. The error model used in the study is one where 
the system errors are assumed to be dominated by point 
measurement noise and data quantization noise. The 
effects of other, more controllable, systematic error 
sources ate not considered in this analysis. 


The superiority of the gravity gradient measuring system 
over the Doppler tracking system found in the studies has been 
corroborated by concurrent studies done elsewhere, such as the 
recent work by Ananda, Lorell and Flury (4] who also compared 
the two systems using a different simulation procedure and got 
essentially the same results. 

Other recent studies also raise some warning about possible 
problems with satellite-to-satellite Doppler tracking systems. 
A recent NASA/Goddard Contract report by Estes and Lancaster (5] 
on a low-low Earth Gravsat mission found that: 

“Recovery of local sets of density blocks from 
long data arcs was found not to be feasible due to 
strong aliasing effects." 

similar studies are underway by Estes on high-low Earth orbit 
satellite-to-satellite Doppler tracking systems which are closer 
to the planned LPO mission configurations. Preliminary results 
indicate that the aliasing problem is less severe for the high- 
low case (factors of 4 to 10 better) but it is still 'nere. It 
would be important for the high-low aliasing ocudies to be 
repeated for typical LPO mission profiles to determine the real 
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accuracy obtainable from an LPO-Relay Doppler tracking system. 
For it may turn out that the data reduction accuracy, when lim- 
ited by aliasing problems, will not permit a determination of 
the mascon densities to the accuracies that would be implied by 
a naive extrapolation of the +1 mm/sec measurement accuracy of 
the Doppler velocity tracking system. 



TASK 4 - SPACECRAFT/GR AD IOMETER INTERFACE CHARACTERISTICS 
AND CONSTRAINTS 

Our work on Task 4 involved the identification of the 
spacecraf t/gradiometer interface characteristics and con- 

straints. These included the following: 

(A) Weight 

(B) Power 

(C) Mechanical Envelope 

(D) Thermal Control 

(E) Science Data Rates/Format 

(F) Engineering Data Rates/Format 

(G) Command (Rate and Format) 

(H) Mechanical Alignment 

(I) Attitude Control (Rate and Acceleration) 


The details of the interface characteristics and con- 
straints as well as a physical and functional description of an 
RGG system designed for the LPO spacecraft are contained in 
Appendix B - Functional and Physical Description of the Gravity 
Environment Measurement System (GEMS) . The general conclusions 
of the study are that the RGG sensor is compatible with the LPO 
spacecraft designs presently under consideration. The detailed 
conclusions can be found in the Appendix. 
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TASK 5 - INCREASED SCIENCE RETURN VS. NUMBER OF SENSORS, 

SPACECRAFT MANEUVERS AND DATA REDUCTION COMPLEXITY 

Task 5 involved an evaluation of the increase in science 
return as a function of the number of sensors (number of axes of 
gravity data) , possible spacecraft maneuvers, and complexity of 
science data reduction. This was done with a computer simula- 
tion of the output of typical RGG sensors operating in an LPO 
spacecraft undergoing orbital and attitude motions that might be 
expected in a typical mission profile. (The computer Fortran 
program is printed in Appendix C.) 

To keep the study within the bounds consistent with the 
extent of this effort, a number of simplifications to the simu- 
lation were made: 

First, the Moon and the LPO orbit were both assumed to be 
insttially fixed in space. The motion of the Moon about the 
Earth and the Earth about the Sun, as well as the effects of the 
other masses in the Solar system were neglected. 

Second, the simulation assumed that the LPO orbit was 
affected only by the primary mass of the moon, and was not 
affected by the rest of the Solar system masses or by the test 
mascons. Since we were not going to calculate the Doppler velo- 
city changes for comparison with the gravity gradient signals we 
did not include the effect of the small test mascons on the 
orbit since the altitude changes (1-5 m) are negligible for the 
test mascon masses that we used in the simulation. 

Third, we assumed that the LPO orbit inclination was xact- 
ly polar. Although it is known that the LPO inclinatio. will 
probably never be greater than 85 degrees because of orbit 
instabilities, we felt that the errors and data reduction prob- 
lems introduced by the non-polar orbit are not of major concern 
and that a simulation carried out assuming a polar orbit would 
show the effects of the various errors that were of concern. 

Fourth, only a few small test mascons were placed in the 
simulation to check the detection sensitivity and resolution of 
the RGG system. We could have attempted to emplace the complex 
multimass system developed in the Task 3 study into this simula- 
tion, but we felt that those results and conclusions could be 
reasonably applied to this simulation if we could show compara- 
ble sensitivity and resolution results for a few point masses. 

Fifth, the error sources introduced into the simulation 
were assumed to be independent. Orbital eccentricity could be 
changed without affecting attitude offsets, and the attitude 
angular rates and pointing errors were both assumed small enough 
so that they could also be considered independent. 
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However, to keep the simulation realistic and to exercise 
the capabilities of the Gravity Gradient measurement system, a 
number of things were not simplified. 

First, the altitude and eccentricity of the LPO orbit was 
made variable. Since a major concern is the change in the back- 
ground bias of the Lunar field as the LPO spacecraft varies in 
altitude, it was important to examine the effects of this back- 
ground bias shift on the data reduction process as it attempted 
to resolve the test mascons. 

Second, although the Moon was assumed to be inertially 
fixed in space, the simulation has the Moon rotating at its 
sidereal angular rate. This rotational motion moves the mascons 
under and past the LPO orbital ground track and allows a deter- 
mination of the cross track resolution that is obtainable from 
comparison of successive orbital passes. 

Third, rather than calculating the gravity gradient tensor 
in a lunar surface fixed frame, the gradient was calculated at 
the position of the spacecraft in an inertially fixed frame. 
This inertially fixed frame rotates with respect to the verti- 
cally oriented spacecraft reference frame once each orbit, which 
rotates the tensor measured by the gravity gradient sensor twice 
each Lunar orbit. Once calculated in the inertial reference 
frame, the gravi ;y gradient tensor is then rotated by the angu- 
lar rotation of the spacecraft to put it into a Lunar vertical 
reference frame. This insures that any subtle mixing effects of 
the local vertical reference system with the ellipticity of the 
orbit will be made evident. 

Fourth, the simulation had provision for inserting both 
angular rate and pointing errors into all three axes of the 
spacecraft measurement frame. These could either be constant 
values or time varying errors. 

Fifth, the simulation provides a multiplicity of outputs. 
These include calculations of the six outputs of the three orth- 
ogonal gravity gradiometer instruments, the six components of 
the gravity gradient tensor field as calculated from the instru- 
ment outputs, the trace of the gravity gradient tensor (which 
should be twice the square of the total spacecraft angular 
rate) , as well as the spacecraft altitude and elapsed time. 

The results of the simulation are presented in Appendix D - 
Lunar Polar Orbiter Gravity Gradient Experiment Simulation 
Plots, by Dr. Robert L. Forward. 

The conclusions that can be drawn from the simulation are: 
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Spacecraft altitude variations due to an elliptical orbit 
may cause some loss of sensitivity and resolution, but the bias 
shifts introduced should not prevent the extraction of the high- 
er frequency gravity data. 

Spacecraft attitude tilts cause a bias shift, and should be 
rapidly ecognizable and measured to milliradian accuracy by 
their appearance in the cross gradient terms. 

Periodic attitude errors caused by malfunction of the atti- 
tude control system will cause significant data reduction errors 
if the attitude error amplitudes are large (greater than 1 
degree) i.nd of rapid period (less than 5 minutes). (Since they 
are tim< varying, they will show up as angular tilts in the 
cross gr; dient outputs and as angular rate gradients of the same 
periodic.ty but quadrature phase in the trace.) 

A single RGG sensor oriented along the orbital axis to mea- 
sure the gravity gradients in the plane of the orbit can use the 
lunar orbit track to track spacing of about 30 km to obtain good 
resolution of mascons in both the along track and cross track 
directions. The cross gradient output of the RGG should be 
integrated with the spacecraft velocity to obtain a contour map 
of tl e vertical gravity field, and then the higher resolution 
(but slightly distorted) principal gradient difference output of 
the RGG should be used to enhance the resolution of the data set 
for geophysical interpretation purposes. 




TASK 6 - POTENTIAL SENSOR SIMPLIFICATIONS FOR ZERO-G 

Task 6 involved evaluating the existing gravity gradiometer 
design to identify those aspects of the instrument which can be 
modified or simplified to reduce instrument costs due to opera- 
tion of the instrument in a zero (0) g environment. 

The RGG sensor design was first examined for potential mod- 
ifications and simplifications with the realization that the 
following RGG error. sources would contribute significantly less 
noise (or none at all) in the zero-g environment of the Lunar 
Polar Orbiter, relative to their noise contributions in the 
earth's g field environment: 

(a) Differential mass unbalance g sensitivity 

(b) Dynamic mass unbalance g sensitivity 

(c) Anisoelastic g-squared sensitivity 

(d) Skew misalignment transverse vibration sensitivity 

(e) Axial vibration-to-torsion sensitivity 

(f) Sum mode mismatch spin variation sensitivity 


If all of these sensitivities were significantly reduced by 
operation in the LPO, then the special designs, devices and tol- 
erance call-outs that were required to keep these sensitivities 
under control in the earth environment could possibly be relaxed 
for the Lunar orbital environment. Some possible cost-saving 
factors that would result would be: 

(a) Elimination of the fine mass balance tubes and their 
control electronics. 

(b) Omission of the isoelastic cut-outs in the sensor arm 
plates. 

(c) Combination of the inner and outer rotor shells into a 
single shell which supports the arm/end mass/pivot structure as 
well as the spin bearings, drag cups and rotary transformers. 

(d) Relaxation of the high tolerances on the fabrication 
of the spin bearings (eliminating the costly final hand lapping 
operation) . 

(e) Reduction in time required for instrument assembly, 
test and balancing. 
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With these possibilities in mind, the analyst who carried 
out the error analyses for the earth bound RGG was asked to 
reexamine his equations and assumptions for an orbital case. 
The relaxation from one g to zero g was found to eliminate the 
need for the fine mass balance tubes and the special isoelastic 
cut-outs in the arm plates, producing some cost savings. Howev- 
er, the largest remaining error source was the skew axis misa- 
lignment driven by transverse vibrations coupling in through the 
misalignment between the spin axis and the sensor torsional 
axis. The transverse vibrations can come from an external 
source or can be generated internally in the spin bearings. 
Although bearing vibrations are usually driven by the g loading 
on the bearing, they can also come from bearing fabrication err- 
ors. Thus, until further work produces a different conclusion, 
we do not feel it is prudent to recommend any further cost 
reductions that could be obtained from relaxing the fabrication 
tolerances on the spin bearings or the cruciform sensing struc- 
ture, or the removal of the provision for aligning the spin axis 
with the torsional axis (which involves differential adjustment 
of the inner rotor holding the sensor/end mass/pivot structure 
and the outer rotor holding the spin bearings) . 
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Appendix A - Gravity Gradient Mapping from the LPO - Science 
Return Simulation Study 

Appendix B - Functional and Physical Description - 

Gravity Environment Measurement System (GEMS) 

Appendix C - LPOSIM Fortran Program - Simulates the 
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Simulation Plots 






LPO GRAVITY GRADIENT EXPERIMENT FORTRAN PROGRAM 


C 

C GRADZZ . FOR 
C 

C CALCULATES THE VERTICAL GRADIENT OF THE VERTICAL GRAVITY 
C OF A SET OF MASS POINTS WHOSE CALCULATED VERTICAL 
C GRAVITY FIELD CLOSELY SIMULATES THE LUNAR GRAVITY FIELD 
C MEASURED BY APOLLO 16 SUBSATELLITE TRACKING DATA FOR 
C THE FRONTSIDE REGION FROM -50 TO +30 DEGREES 
C LONGITUDE AND +4 TO +9 DEGREES LATITUDE. 

C REFERENCE* 

C W. L. SJOGREN , R. N. WIMBERLY AND W. R. WOLLENHAUPT 

C "LUNAR GRAVITY VIA THE APOLLO 15 AND 16 SUBSATELLITES" 

C THE MOON, 9, 115-128 (1974). 

C 

C 30 JUNE 1976 VERSION 
C 

C ORIGINAL PROGRAM WRITTEN BY MR. ROGER W. LOWE 
C HUGHES AIRCRAFT COMPANY 
C 

C REFER QUESTIONS OR COMMENTS TO: 

C DR. ROBERT L. FORWARD 

C HUGHES RESEARCH LABORATORIES 

C 3011 MALIBU CANYON ROAD 

C MALIBU, CA 90265 

C (213)456-6411 X280 

C 

C ONE SET OF INPUT DATA ARE THE MASCONS STORED IN 
C MASCON.DAT THAT ARE IN A 15X181 ARRAY- WITH 15 KM SPACING 
C ON THE LUNAR SURFACE. 

C THE PROGRAM USES A 9X9 SUBARRAY OF MASCONS TO CALCULATE 
C THE GRADIENT AT EACH POINT. 

C 

C THE OTHER SET OF INPUT DATA ARE THE LARGE DOMINANT MASCONS 
C WHOSE MASS AND POSITION ARE STORED IN MASXYZ.DAT. 

C 

C THE UNITS USED ARE: 

C DISTANCES IN KILOMETERS 

C MASSES IN MASCONS (10**14 KG) 

C GRAVITY GRADIENT IN EOTVOS (10**-9 SEC**-2) 

C 

C IN THESE UNITS T8E NEWTONIAN GRAVITATIONAL CONSTANT IS 
C G-6670 EOTVOS KM**3/MASCON 


preceding page blank not 
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DIMENSION THE DATA AND OPEN THE FILES 

DIMENSION RM(201,35) ,GRDZZ (15) 

DIMENSION IX ( 30) ,IY(30) ,A(30) ,W(30) 

OPEN (UNIT-20, FILE=' MASCON.DAT' ) 

OPEN (UNIT=21 , FILF= ’ MASXYZ . DAT ' ) 

OPEN (UN IT=2 2 , F I LL= 'GRADZZ . DAT ' ) 

READ IN THE SMALL MASCON GRID POINTS 

DO 10 1=11,191 

READ (20,1) (RM(I,J) ,J=11,25) 

1 FORMAT (15F5.1) 

10 CONTINUE 

FOLD OUT A COPY OF THE OUTER 10 POINTS OF THE MASCON 
GRID TO GENERATE A LARGER GRID FOR THE CALCULATION WITH 
A SMOOTH TRANSITION AT THE MEASUREMENT BOUNDARY 


DO 20 1=11,191 

DO 15 J = 1 ,10 

RM ( I , J ) =RM ( 1 , 21-J) 

15 CONTINUE 

DO 20 J=26 , 35 
RM ( I , J) =RM ( I , 51-J) 

20 CONTINUE 

DO 40 J = ll , 25 

DO 30 1=1,10 

RM ( I , J ) = RM ( 2 1 - 1 , J ) 

30 CONTINUE 

DO 40 1=192,201 
RM(I,J)=RM(383-I,J) 

40 CONTINUE 

DO 50 1*1,10 
DO 50 J=1,10 
RM ( I , J) =RM ( I , 21-J) 

RM ( I , J+25) *RM ( I ,26-J) 

50 CONTINUE j ! 

DO 60 1 = 192,201 ! | ! 

DO 60 J»l,10 I 

RM ( I , J) =RM ( I , 21-J) ! 

RM ( I , J+25) *RM ( I , 26-J) j 

60 CONTINUE 

0 • 
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READ IN THE LARGE MASCON DATA SET 

IX (I) * X POSITION OF MASCON IN UNITS OF 15 KM 
IY (I) * Y POSITION OF MASCON IN UNITS OF 15 KM 
A (I) = DEPTH OF MASCON BELOW SURFACE IN KM 
W(I) = WEIGHT OF MASCON IN UNITS OF 10**14 KG 

DO 70 1=1,30 

READ (21 , 3 , END=80) IX(I) ,IY(I) ,A(I) ,W(I) 

3 FORMAT (1X,I3,1X,I3,1X,F4.0,F7.0) 

70 CONTINUE 
80 NUM-I-1 

SET THE ALTITUDE ABOVE THE LUNAR SURFACE 
H= ALTITUDE IN KILOMETERS 

H=30. 

CALCULATE THE GRAVITY GRADIENT FROM THE 81 SMALLER 
MASCONS IN A 9X9 SUBARRAY ABOUT THE MEASUREMENT POINT. 
X,Y,Z= X, Y, Z DISTANCE FROM MEASUREMENT POINT IN KM 

Z*0.-H 

DO 120 1=11,191 
DO 110 J*ll ,25 
G=0 

DO 100 M— 4,4 
X=M*15. 

DO 90 N»-4,4 
Y=N*15. 

P2«X*X+Y*Y 
R*SQRT ( P2+ Z * Z ) 

G«G+6670.*RM(I+M,J+N)*(2-3*P2/(R**2))/(R**3) 

90 CONTINUE 
100 CONTINUE 

CALCULATE THE GRAVITY GRADIENT FROM THE 
LARGER MASCON SET. 

DO 105 K=1 ,NUM 
X= (IX(K) — ( I— 10 ) ) *15. 

Y=(IY(K)-(J-10))*15. 

Z=-(A(K)+H) 
t >X*X+Y*Y 
R=SQRT(P2+Z*Z) 

G=G+6670.*W(K)*(2-3*P2/(R**2))/(R**3) 

105 CONTINUE 
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STORE THE TOT A I. INTEGRATED GRAVITY GRADIENT 

FOR THAT MEASUREMENT POINT AMD MOVE ON TO THE NEXT 

MEASUREMENT POINT. 

GRDZZ (J-10) =G 
110 CONTINUE 

WRITE (22,4) GRDZZ 
4 FORMAT ( 1 5F5 . 1 ) 

120 CONTINUE 

STOP WHEN FINISHED 

STOP 

END 
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LARGE MASCON DATA SET (MASXYZ.DAT) 


1 

{ 

V 

nL 


This is a set of 20 large mass points at various positions 
and depths below the lunar surface. The first two numbers are 
the x and y positions of the masses in integer units of 15 Km, 
the third number is the depth in Km, and the fourth is the mass 
in units cf 10*14 Kg. 


1 1 


r ; 

i 


li 

0 

D 

D 

D 

D 

0 

1 


003 

013 

050. 

440. 

007 

-02 

065. 

-620. 

016 

014 

045. 

-160. 

018 

-01 

085. 

970. 

027 

014 

45. 

-230. 

033 

009 

45. 

-100. 

043 

Oil 

75. 

440. 

050 

014 

40. 

-200. 

053 

006 

35. 

170. 

060 

003 

60. 

-260. 

062 

016 

30. 

-250. 

067 

Oil 

50. 

800. 

078 

010 

25. 

160. 

084 

015 

55. 

900. 

090 

006 

60. 

460. 

096 

013 

45. 

-470. 

101 

008 

40. 

-420. 

104 

-01 

60. 

540. 

110 

007 

50. 

470. 

117 

010 

75. 

540. 

119 

002 

50. 

-340. 

126 

002 

75. 

-890. 

128 

012 

75. 

-1010. 

148 

007 

60. 

1100. 

160 

007 

75. 

-680. 

166 

-Qi 

85. 

-820. 
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SMALL MASCON DATA SET (MASCON.DAT) 

This is a compacted version of a 15x181 set of mass points 
that are spaced at 15 Km intervals on the Lunar surface. The 
units of the mascon masses are in 10*14 Kg. 


-0.4 -0.5-0. 6-0. 7-0. 7-0. 7-0. 7-0 
-0.2-0. 5-0. 6-0. 8-0. 8-0. 8-0. 9-1 
-0.2-0. 6-1. 2-1. 0-0. 8-0. 8-0. 9-1 
+0.1-0. 6-1. 2-1. -0.8-0. 9-0. 8-0 
+0.5-0.5-0.7-1.0-0.9-0.8-0.7-0 
40.4-0.5-1.0-1.0-0.6-0.7-0.4+0 
+0.3-1. 0-1. 5-0. 6-0. 5-0. 3-0. 2+0 
+0.3-1. 0-1. 1-0. 5-0. 4-0. 3-0. 1-0 
-0 . 4-0 . 9-0 . 4- 0 . 3-0 . 2-0 . 2-0 . 1-0 
-0.7-0. 5-0. 3-0. 1+0. 0-0. 0-0. 2+0 
-0.5-0. 4-0. 1+0. 1+0. 3+0. 1-0. 4-0 
-0 . 5-0 . 2-0 . 0+0 . 6+0 . 8+0 . 8+0 . 5+0 
-0.3-0. 2+0. 5+1. 3+1. 6+1. 5+1. 5+2 
-0.7+0. 4+1. 1+1. 7+1. 7+1. 9+2. 0+1 
-0.2+0. 8+1. 6+1. 9+1. 8+1. 7+1. 6+1 
-0.0+1. 1+2. 0+1. 7+1. 5+1. 7+1. 2-0 
-0.0+0. 7+1. 5+1. 3+1. 6+1. 8+0. 9-0 
+0.2+0. 6+1. 6+1. 0+1. 2+1. 4-1. 0+0 
+0.1+0. >+0.4+0. 3+0. 1-0. 8-1. 3-0 
+0.3+0. 4+0. 1+0. 3-0. 7-2. 1-1. 2-0 
+0.1-0.1-0.3-0.6-2.5-1.2-0.9-0 
+0.2-0.3-0.6-1.0-1.8-1.1-0.9-0 
+0 . 0-0 . 5-0 . 6-1 . 6-1 . 0-0 . 9-0 . 8-0 
+0.0-0.5-0.7-1.1-1.0-0.8-0.8-0 
-0 . 3-0 . 5-0 . 8-0 . 7-0 . 7-0 . 6-0 . 5-0 
+0.0-0.3-0.7-0.8-0.6-0.5-0.4-0 
+0.2-0.3-0.4-0.8-0.6-0.4-0.3-0 
+0.3-0.1-0.3-0.6-0.4-0.3-0.2-0 
+0.7+0. 1-0. 4-0. 4-0. 3-0. 2-0. 1-0 
+0.6+0. 1-0. 4-0. 4-0. 2-0. 1+0. 0+0 
+0 .6+0 . 3-0 . 5-0 . 3-0 . 3-0 . 1+0 .1+0 
+1.0+0. 5-0. 1-0. 3-0. 2-0. 2+0. 0+0 
+1.6+0. 7+0. 2-0. 3-0. 2-0. 2-0. 3-0 
+1.6+3 .0+0. 2+0. i-0. 1-0. 3-0. 7-0 
+2 . 0+1 . 3+1 . 2+0 . 7+0 . 8-0 . 0-0 . 2-0 
+2. 1+1. 6+1. 4+1. 0+0. 7+0. 4-0. 0-0 
+1.9+1. 8+1. 5+1. 3+1. 0+0. 6+0. 1+0 
+2. 1+1. 9+1. 8+1. 6+1. 4+0. 6+0. 6+0 
+2. 3+2. 0+1. 7+1. 6+1. 6+1. 9+1. 6+0 
+2. 3+2. 2+1. 9+2. 2+2. 2+1. 9+1. 7+1 


.7-0. 4+0. 0+1. 1+1. 1+0. 7+0. 4+ 1.6 
.1-1. 0-0. 2-0. 0-0. 0+0. 6-0. 5+0. 5 
.0-0. 4-0. 2-0. 1-1. 4+0. 6-0. 3+0.1 
. 9+0 . 1+0 . 6-0 . 8-0 . 9-0 . 2+0 .4+0.1 
.5+0. 1+0. 6-0. 1-0. 1-0. 6-0. 3+0.0 
. 1+0 . 3+0 . 5+0 . 1+0 . 3+0 . 0+0 .0+0.2 
.3+0. 4+0. 6+0. 6+0. 4+0. 2-0. 1-0. 4 
.1-0. 0+0. 4+0. 6+0. 7+0. 9-0. 4-0. 5 
. 1 -0 . 1+0 . 0+0 . 3-0 . 1-0 . 5-0 . 7-0 . 8 
. 0-0 . 1-0 . 2-0 . 5-1 . 1-1 . 3-1 . 6-1 . 6 
.2+0. 0-0. 3-1. 0-1. 6-1. 3-1. 5-1. 9 
.4+0. 1-0. 4-0. 9-1. 2-1. 3-1. 3-1. 6 
.5+0. 4-0. 5-0. 8-1. 1-1. 1-1. 0-1.0 
.6+0. 5-1. 1-0. 7-0. 8-1. 0-0. 7-0. 7 
.9+0. 6-0. 7-0. 6-0. 6-0. 9-0. 5+0.0 
. 6+0 . 2-0 . 4-0 . 7-0 . 5-0 . 5-0 .2+0.1 
.3+0. 2-0. 4-0. 6-0. 5-0. 4-0. 1+0.1 
.1+0. 5+0. 0-0. 8-0. 5-0. 4-0. 5-0. 4 
.4+0. 0+1. 4-1. 0-0. 8-0. 5-0. 6-0. 6 
. 3-0 . 1+0 . 5-0 . 3-0 . 3-0 . 6-0 . 8-1 . 0 
.4+0 . 0+0 . 3+0 .4-0 . 2-0 . 4-0 . 8-1 . 2 
.4-0. 1+0. 4+1. 2+0. 3-0. 3-0. 7-1.1 
. 4+0 . 0+0 . 5+1 . 0+0 . 5-0 . 2-0 .7-0.8 
. 4-0 . 2-0 . 2+0 . 1+0 . 1 *0 , 1-0 .2-0.3 
.4-0. 5+0. 0+0. 1+0. 1+0. 1+0. 0+0. 2 
. 3-0 . 1+0 . 1+0 . 2+0 . 6+0 . 1-1 . 3+0 . 6 
.1-0. 0+0. 4+0. 5+1. 2+0. 2-0. 6+1. 4 
. 1+0 . 1+0 . 4+0 . 8+1 . 1-0 . 0-1 . 5+1 . 2 
.0+0. 1+0. 4+1. 0+0. 4-0. 1-0. 1+0. 2 
. 1+0 . 4+0 . 4+0 . 6+0 . 3-0 . 3-0 .3-0.1 
. 3+0 . 7+0 . 7+0 . 5-0 . 2-0 . 2-0 . 3-0 . 3 
. 6+1 . 1+0 . 8+0 . 2-0 . 1-0 . 3-0 . 4-0 . 4 
.2-0. 9-0. 3-0. 9-0. 4-0. 5-0. 8-0. 7 
. 6-1 . 2-0 . 6-0 . 2-0 . 6-1 . 1-0 .9-0.8 
. 6-1 . 2-0 . 5-0 . 2-0 . 8-0 . 9-1 . 4+0 . 6 
.4-0. 6-0. 2-0. 3-0. 6-1. 2-0. 5+0. 4 
.1+0. 1+0. 6-0. 2-0. 6-0. 8-0. 9+1. 7 
. 3+0 . 2+0 . 1-0 . 1-0 . 1-0 . 3-0 .0+0.7 
. 8+0 . 3-0 . 1+0 . 0-0 . 1-0 . 0+0 .5+0.7 
. 2+0 . 3+0 . 1-0 . 2-0 . 2-0 . 5+0 . 6+0 . 5 
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+ 2 . 7 + 2 . 1 + 1 . 7 + 1 . 1 + 1 . 3 + 1 . 3 + 1 . 0 + 0 . 5 - 0 . 7 + 0 . 1 + 0 . 1 + 0 . 1 + 0 . 5 + 1 . 0 + 0. 4 
+ 2 . 8 + 2 . 2 + 1 . 5 + 1 . 2 + 1 . 1 + 1 . 4 + 0 . 9 + 0 . 4 - 0 . 8 - 0 . 3 - 0 . 1 + 0 . 1 + 0 . 5 + 0 . 6 + 0. 5 
+ 3 . 0 + 1 . 8 +-. 3 + 0 . 3 + 0 . 6 + 0 . 7 + 0 . 7 + 0 . 3 - 0 . 7 - 0 . 3 + 0 . 1 + 0 . 2 + 0 . 5 + 0 . 2 + 0. 2 
+ 2 . 7 + 1 . 5 + 0 . 5 + 0 . 1 + 0 . 1 + 0 . 3 + 0 . 3 - 0 . 3 - 0 . 3 + 0 . 1 + 0 . 1 + 0 . 1 + 0 . 1 + 0 . 3 - 0.1 
+1 . 9+1 . 1+0 . 5+0 . 1-0 . 0+0 . 2+0 . 3+0 . 2+0 . 3+0 . b +0 . 1-0 . 0+0 . 9-0 . 2-1 . 0 
+ 1 . 7 + 0 . 6 + 0 . 6 + 0 . 2 + 0 . 1 + 0 . 2 + 0 . 3 + 0 . 3 + 0 . 5 + 0 . 9 + 0 . 3 - 0 . 1 - 0 . 5 - 0 . 7 - 2.1 
+1 . 0+0 . 8+0 . 5+0 . 2+0 . 1+0 . 2+0 . 3+0 . 4+0 . 8+1 . 3+0 . 1-0 . 4-0 . 7-2 . 3-2 . 0 
+ 0 . 8 + 0 . 8 + 0 . 8 + 0 . 8 + 0 . 4 + 0 . 2 + 0 . 3 + 0 . 4 + 1 . 1 + 0 . 4 - 1 . 2 - 0 . 8 - 1 . 2 - 0 . 8 - 1. 3 
+ 1 . 0 + 1 . 0 + 1 . 0 + 0 . 7 + 0 . 4 - 0 . 1 - 0 . 1 + 0 . 3 + 0 . 5 + 0 . 6 - 0 . 1 - 0 . 9 + 0 . 6 + 0 . 3 - 0. 2 
+ 1 . 1 + 1 . 3 + 1 . 4 + 0 . 9 + 0 . 4 - 0 . 1 + 0 . 2 + 0 . 4 + 1 . 0 + 0 . 5 - 0 . 1 - 0 . 6 + 0 . 4 - 0 . 0 + 0. 5 
+ 1 . 1 + 1 . 2 + 1 . 2 + 1 . 1 + 0 . 2 - 0 . 7 + 2 . 2 + 1 . 2 + 1 . 3 + 0 . 5 + 0 . 7 - 1 . 0 - 0 . 2 - 0 . 3 + 0. 6 
+ 1 . 2 + 1 . 3 + 1 . 0 + 0 . 6 - 0 . 0 - 0 . 2 + 0 . 3 + 0 . 5 + 1 . 1 + 1 . 8 + 0 . 2 - 0 . 8 - 0 . 5 - 0 . 7 - 0.3 
+ 1 . 2 + 1 . 2 + 1. 1 ‘- 0 . 5 - 0 . 3 - 0 . 8 - 0 . 6 + 0 . 2 + 1 . 0 + 1 . 4 + 0 . 3 - 0 . 7 - 1 . 9 - 1 . 1 - 1.0 
+ 1 . 0 + 0 . 6 + 0 . 7 + 1 . 4 + 0 . 7 - 0 . 7 - 0 . 4 - 1 . 4 + 0 . 3 + 1 . 2 + 0 . 7 + 1 . 1 - 1 . 0 - 2 . 0 - 0.7 
+ 1 . 0 - 0 . 3 - 1 . 6 + 0 . 7 + 0 . 9 + 0 . 0 - 0 . 3 - 0 . 4 - 0 . 6 + 0 . 3 + 0 . 7 + 0 . 6 + 0 . 0 - 1 . 0 + 0.3 
- 0 . 3 - 0 . 9 - 1 . 4 - 1 . 0 - 1 . 0 + 1 . 1 - 0 . 2 - 0 . 3 - 0 . 1 + 0 . 3 + 0 . 8 + 1 . 3 + 0 . 3 - 0 . 0 + 1.1 
- 0 . 3 - 0 . 4 - 1 . 7 - 0 . 7 - 0 . 6 - 1 . 3 - 0 . 4 - 0 . 3 + 0 . 1 + 0 . 4 + 0 . 9 + 0 . 8 + 0 . 5 + 0 . 2 rl.O 
- 0 . 3 - 0 . 2 - 0 . 6 - 1 . 0 - 0 . 6 - 0 . 4 - 0 . 3 - 0 . 1 + 0 . 1 + 0 . 3 + 1 . 1 + 0 . P + 0 . 7 + 0 . 7 + 0.8 
+0 . 0+0 . 9-0 . 4-0 . 4-0 . 5-0 . 3-0 . 2-0 . 1+0 . 1+0 . 2+0 . 7+0 . i +0 . 8+0 . 6+0 . 7 
- 1 . 4 - 0 . 5 + 0 . 0 - 0 . 4 - 0 . 5 - 0 . 5 - 0 . 3 - 0 . 1 + 0 . 2 + 0 . 4 + 1 . 0 + 0 . 9 + 1 . 0 + 1 . 0 + 0. 5 
- 0 . 3 + 0 . 2 + 0 . 9 - 0 . 4 - 0 . 8 - 0 . 8 - 0 . 7 - 0 . 1 - 0 . 4 + 0 . 2 + 2 . 6 - 0 . 4 + 1 . 4 + 1 . 4 + 0. 8 
- 0 . 1 + 0 . 1 + 0 . 0 - 0 . 5 - 1 . 5 - 1 . 4 - 0 . 9 - 1 . 2 - 2 . 6 + 4 . 6 + 2 . 6 - 0 . 0 + 1 . 8 + 2 . 2 - 0. 8 
+ 0 . 3 - 0 . 2 - 0 . 4 - 0 . 8 - 1 . 0 - 1 . 4 - 1 . 8 - 0 . 7 + 1 . 8 + 3 . 0 + 4 . 3 + 0 . 8 + 3 . 9 - 0 . 6 - 3. 9 
+ 0 . 9 + 0 . 2 - 0 . 7 - 1 . 2 - 1 . 0 - 1 . 2 + 1 . 9 - 0 . 4 - 1 . 0 + 2 . 0 + 2 . 4 + 2 . 4 + 4 . 2 - 6 . 6 - 1. 8 
+ 0 . 3 + 0 . 4 - 1 . 1 - 1 . 5 - 1 . 4 - 2 . 0 + 1 . 7 - 2 . 9 + 2 . 1 + 1 . 0 + 0 . 3 + 1 . 9 + 3 . 9 - 0 . 2 - 1. 6 
+ 0 . 2 + 0 . 0 - 0 . 5 - 1 . 4 - 1 . 8 - 0 . 8 + 0 . 2 + 0 . 7 + 3 . 4 + 0 . 4 - 0 . 6 - 1 . 2 + 3 . 7 + 2 . 0 - 2.0 
+ 0 . 1 - 0 . 1 - 0 . 6 - 1 . 7 - 2 . 3 - 1 . 4 + 0 . 2 + 1 . 3 + 1 . 6 - 0 . 2 - 1 . 5 - 0 . 0 + 1 . 3 + 3 . 7 + 0. 9 
+ 0 . 0 - 0 . 2 - 0 . 7 - 1 . 6 - 2 . 6 - 1 . 9 + 0 . 7 + 2 . 0 + 1 . 4 + 0 . 3 + 0 . 1 + 0 . 5 + 1 . 7 + 3 . 6 + 1. 7 
+ 0 . 0 - 0 . 2 - 0 . 7 - 1 . 5 - 2 . 5 - 2 . 3 + 0 . 2 + 1 . 9 + 1 . 0 + 1 . 1 + 1 . 0 + 2 . 0 + 1 . 9 + 2 . 4 - 0. 4 
+ 0 . 2 - 0 . 1 - 0 . 7 - 1 . 4 - 2 . 3 - 1 . 4 - 0 . 6 + 0 . 7 + 0 . 6 + 1 . 1 + 0 . 9 + 2 . 3 - 0 . 9 + 1 . 2 - 0. 8 
+ 0 . 0 - 0 . 2 - 0 . 6 - 1 . 5 - 3 . 0 - 0 . 1 - 2 . 5 + 1 . 5 + 0 . 6 + 1 . 8 + 4 . 9 + 0 . 9 + 0 . 7 - 0 . 4 - 0.8 
+ 0 . 0 - 0 . 1 - 0 . 2 - 0 . 5 - 0 . 7 - 1 . 1 - 1 . 4 - 0 . 6 - 0 . 5 - 1 . 9 + 1 . 0 - 0 . 4 - 0 . 5 - 0 . 3 - 0. 5 
+ 0 . 0 - 0 . 3 + 0 . 1 + 0 . 3 + 0 . 8 - 1 . 3 - 1 . 3 - 0 . 6 - 0 . 6 - 1 . 1 + 2 . 4 - 0 . 2 + 0 . 0 - 0 . 3 + 0.0 
- 0 . 2 - 0 . 5 + 0 . 5 + 1 . 4 + 1 . 7 + 2 . 8 - 1 . 6 - 0 . 6 - 0 . 0 + 2 . 8 + 1 . 6 + 0 . 3 - 0 . 2 - 0 . 4 + 0.1 
- 0 . 0 + 0 . 4 + 1 . 0 + 1 . 7 + 2 . 4 + 4 . 9 - 0 . 0 - 1 . 3 - 0 . 1 + 4 . 3 + 2 . 7 - 0 . 7 - 0 . 2 - 0 . 2 + 0. 6 
- 0 . 0 + 0 . 8 + 1 . 7 + 2 . 2 + 2 . 0 + 2 . 9 + 1 . 8 + 3 . 7 + 2 . 9 + 2 . 6 + 2 . 7 - 0 . 1 - 0 . 9 + 0 . 1 + 1.1 
- 0 . 3 + 2 . 3 + 1 . 8 + 1 . 6 + 1 . 4 + 2 . 5 + 2 . 2 + 2 . 5 + 0 . 8 - 1 . 0 + 1 . 1 - 0 . 2 - 0 . 3 - 0 . 0 + 1. 8 
- 1 . 5 - 0 . 7 + 0 . 7 + 1 . 4 + 1 . 3 + 1 . 2 + 1 . 5 + 0 . 7 - 1 . 6 - 1 . 0 - 0 . 5 + 0 . 3 - 0 . 1 + 0 . 7 + 2.1 
- 1 . 0 - 1 . 0 + 0 . 1 + 0 . 5 + 0 . 8 + 0 . 8 + 0 . 3 - 1 . 3 + 0 . 7 - 1 . 3 + 0 . 4 + 0 . 9 + 0 . 4 + 0 . 4 + 2. 7 
- 1 . 1 - 1 . 2 - 0 . 1 + 0 . 1 + 0 . 3 + 0 . 3 + 0 . 1 - 0 . 8 - 1 . 4 + 1 . 4 + 0 . 9 + 1 . 8 + 1 . 4 + 2 . 0 + 3. 6 
- 1 . 0 - 1 . 5 - 0 . 2 + 0 . 0 - 0 . 1 + 0 . 1 + 0 . 6 - 1 . 1 - 0 . 8 + 0 . 2 + 1 . 2 + 0 . 9 + 1 . 0 + 1 . 6 + 2. 8 
- 1 . 1 - 1 . 6 + 0 . 2 - 0 . 1 + 0 . 3 + 0 . 7 + 0 . 3 - 1 . 0 - 1 . 3 - 0 . 1 + 0 . 9 + 0 . 4 - 0 . 2 - 0 . 8 + 1.2 
- 1 . 2 - 1 . 9 - 0 . 3 - 0 . 2 + 1 . 2 + 0 . 9 + 0 . 1 - 0 . 8 - 1 . 8 - 0 . 2 - 0 . 1 - 0 . 0 - 0 . 0 + 0 . 2 + 1. 2 
- 1 . 7 - 0 . 5 - 0 . 1 + 0 . 8 + 1 . 0 + 0 . 5 + 0 . 2 - 0 . 5 - 1 . 4 - 1 . 5 - 0 . 6 - 0 . 2 + 0 . 1 - 0 . 7 - 0.0 
- 0 . 7 - 0 . 2 + 0 . 6 + 0 . 9 + 0 . 5 + 0 . 2 + 0 . 1 - 0 . 5 - 1 . 6 - 1 . 6 - 0 . 9 - 0 . 4 + 0 . 5 + 0 . 7 + 0. 3 
- 0 . 4 + 0 . 3 + 1 . 1 + 0 . 2 - 0 . 0 - 0 . 1 - 0 . 3 - 0 . 5 - 1 . 1 - 1 . 7 - 1 . 5 - 0 . 9 + 0 . 5 + 0 . 4 - 0. 2 
+ 0 . 1 + 0 . 1 + 0 . 3 - 0 . 1 - 0 . 2 - 0 . 4 - 0 . 5 - 0 . 7 - 1 . 0 - 1 . 7 - 2 . 6 + 0 . 2 - 0 . 8 + 1 . 3 + 1.0 
+ 0 . 1 + 0 . 2 + 0 . 0 - 0 . 5 - 0 . 9 - 0 . 5 - 0 . 7 - 0 . 9 - 0 . 3 - 0 . 1 - 3 . 5 - 1 . 7 + 1 . 5 + 1 . 7 + 2. 9 
+ 0 . 1 + 0 . 2 - 0 . 1 - 1 . 2 - 1 . 3 - 0 . 6 - 0 . 8 - 0 . 9 - 0 . 8 - 2 . 8 - 3 . 8 - 2 . 4 - 1 . 7 + 0 . 6 + 2.0 
+0.3+0. 1-0. 4-2. 5+3. 7-0. 7-0. 8-1. 1-1. 1-1. 3-3. 3-1. 7-1. 1-1. 8+1. 3 
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+ 0 . 8 + 1 . 5 + 1 . 7 + 1 . 6 + 0 . 7 + 0 . 3 - 0 . 1-0 
+ 2 . 1 + 2 . 5 + 2 . 2 + 1 . 9 + 1 . 1 + 0 . 7 + 0 . 1-0 
+ 2 . 5 + 2 . 8 + 2 . 3 + 2 . 2 + 1 . 8 + 1 . 4 + 0 . 1-0 
+ 2 . 6 + 2 . 6 + 2 . 6 + 2 . 3 + 1 . 9 + 0 . 8 - 1 . 0-1 
+ 2 . 5 + 2 . 3 + 2 . 3 + 2 . 0 + 1 . 5 + 0 . 5 - 1 . 6-1 
+ 1 . 8 + 2 . 1 + 2 . 3 + 3 . 1 + 1 . 1 + 0 . 3 + 2 . 0-1 
+ 0 . 5 + 1 . 7 + 1 . 6 + 1 . 3 + 0 . 7 - 1 . 8 - 0 . 4-0 
- 0 . 0 + 2 . 1 + 1 . 0 - 0 . 8 + 1 . 8 - 2 . 1 + 1 . 3-0 
- 0 . 6 + 1 . 4 + 0 . 5 - 0 . 7 + 1 . 5 - 0 . 3 - 1 . 6-1 
- 0 . 0 - 0 . 2 + 0 . 1 + 0 . 1 - 0 . 0 - 0 . 6 - 0 . 7-0 
- 0 . 1 - 0 . 1 - 0 . 7 - 0 . 0 + 0 . 5 + 0 . 8 + 0 . 4+0 
- 0 . 1 - 0 . 2 - 0 . 2 - 0 . 1 + 1 . 3 + 1 . 0 + 0 . 4-0 
- 0 . 2 - 0 . 3 - 0 . 3 - 0 . 6 - 0 . 1 - 0 . 1 - 0 . 5+1 
- 0 . 3 - 0 . 8 - 1 . 5 - 0 . 3 + 0 . 7 + 0 . 5 + 0 . 3+0 
+ 0 . 2 - 0 . 8 - 0 . 7 - 0 . 3 - 0 . 1 + 0 . 2 + 0 . 1-0 
+ 0 . 3 - 1 . 4 - 1 . 0 + 0 . 5 + 0 . 0 + 0 . 0 + 0 . 6-0 
- 0 . 8 - 1 . 8 - 1 . 6 + 0 . 1 + 0 . 2 - 1 . 3 + 0 . 3+0 
- 1 . 1 - 2 . 5 - 1 . 3 - 0 . 7 - 1 . 1 - 0 . 2 + 0 . 2+0 
- 1 . 6 - 2 . 0 - 1 . 8 - 0 . 4 - 0 . 1 + 0 . 1 + 1 . 1+0 
- 0 . 4 - 2 . 4 - 1 . 3 - 0 . 5 + 0 . 1 - 0 . 3 + 0 . 4+0 
- 1 . 9 - 1 . 3 - 0 . 8 - 0 . 5 + 0 . 1 - 0 . 2 + 0 . 0+0 
- 1 . 2 - 0 . 8 - 0 . 6 - 0 . 1 - 0 . 3 - 0 . 3 - 0 . 2-0 
- 0 . 7 - 0 . 5 - 0 . 4 - 0 . 3 - 0 . 3 - 1 . 4 - 0 . 6-0 
- 0 . 4 - 0 . 2 - 0 . 3 - 0 . 3 - 0 . 2 + 0 . 4 - 0 . 3-0 
- 0 . 2 + 0 . 0 - 0 . 2 - 0 . 2 + 0 . 0 + 0 . 0 - 0 . 3-0 
+ 0 . 0 + 0 . 0 + 0 . 0 + 0 . 0 + 0 . 2 - 0 . 2 - 0 . 4-0 
- 0 . 0 + 0 . 0 - 0 . 2 - 0 . 4 - 0 . 3 - 0 . 1 - 0 . 4-0 
+ 0 . 5 + 0 . 1 - 0 . 3 - 0 . 5 - 0 . 3 + 0 . 5 - 0 . 3-0 
- 0 . 2 - 0 . 5 - 0 . 6 - 0 . 8 - 0 . 3 + 0 . 0 - 0 . 2-0 
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- 1 . 1 - 1 . 2 - 1 . 2 - 2 . 1 - 0 . 3 - 0 . 3 - 0 . 6-0 
- 1 . 3 - 1 . 5 - 1 . 6 - 1 . 1 - 0 . 8 - 0 . 7 - 0 . 7-0 
- 1 . 8 - 1 . 8 - 1 . 4 - 0 . 8 - 0 . 8 - 0 . 7 - 0 . 7-1 
- 2 . 0 - 1 . 9 - 1 . 3 - 0 . 9 + 0 . 2 - 0 . 7 - 0 . 7-0 
- 2 . 1 - 2 . 2 - 1 . 6 - 0 . 5 - 0 . 1 - 0 . 4 - 0 . 5-0 
- 2 . 1 - 2 . 4 - 2 . 2 - 0 . 6 - 0 . 3 - 0 . 5 - 0 . 5-0 
- 2 . 4 - 2 . 7 - 2 . 5 + 0 . 0 - 0 . 7 - 0 . 6 - 0 . 6-0 
- 2 . 5 - 2 . 8 - 1 . 8 - 0 . 2 - 0 . 7 - 0 . 6 - 0 . 8-0 
- 2 . 3 - 2 . 3 - 1 . 6 - 0 . 4 - 0 . 7 - 0 . 6 - 0 . 5-0 
- 2 . 3 - 1 . 9 - 1 . 3 - 0 . 5 - 0 . 7 - 0 . 7 - 0 . 9-0 
- 2 . 0 - 1 . 7 - 0 . 8 - 0 . 7 - 0 . 8 - 0 . 6 - 0 . 6+0 


. 2 - 0 . 3 - 0 . 5 - 0 . 5 - 0 . 6 - 0 . 6 - 0 . 7 - 0.4 
. 9 - 0 . 3 - 0 . 3 - 0 . 1 + 0 . 0 + 0 . 2 - 0 . 8 - 0.9 
. 1 - 0 . 2 - 1 . 0 - 0 . 1 + 0 . 1 + 0 . 4 + 1 . 1 + 0 . 3 
. 0 + 0 . 2 + 0 . 9 + 0 . 1 + 0 . 6 + 0 . 9 + 0 . 9 + 1.0 
. 8 - 1 . 6 - 0 . 4 + 0. 3 + 1 . 5 + 1 . 3 + 1 . 2 + 1. 3 
. 1 - 1 . 1 - 0 . 8 + 0 . 3 + 1 . 1 + 1 . 4 + 1 . 7 + 1. 6 
. 3 - 0 . 6 - 0 . 8 + 0 . 4 + 1 . 1 + 1 . 7 + 2 . 2 + 2. 3 
. 4 - 0 . 4 - 0 . 9 + 0 . 5 + 1 . 5 + 1 . 6 + 2. 1 + 2. 3 
. 7 + 0 . 6 - 0 . 4 + 0 . 6 + 1 . 1 + 1 . 7 + 1 . 9 + 2. 4 
. 4 + 0 . 1 + 0 . 2 + 0 . 5 + 1 . 0 + 1 . 6 + 1. 8 + 1. 9 
. 5 + 0 . 5 + 0 . 3 + 0 . 4 + 0 . 7 + 0. 7 + 1 . 3 + 1 .8 
. 2 + 0 . 1 + 0 . 6 + 0 . 7 + 0 . 4 + 0. 7 + 1. 0 + 1. 4 
. 0 + 0 . 7 + 0 . 0 - 0 . 1 + 0 . 1 + 0 . 3 + 0. 7 + 1. 2 
. 1 - 0 . 5 - 0 . 6 - 1 . 0 - 0 . 8 - 0 . 7 - 0 . 4 + 0. 7 
. 6 - 0 . 8 - 1 . 7 - 1 . 9 - 0 . 9 - 0 . 4 - 0 . 2 - 0.0 
. 3 - 0 . 3 - 1 . 1 - 1 . 1 - 0 . 8 + 0 . 1 - 0 . 1 + 0.0 
. 4 - 0 . 5 - 1 . 5 - 2 . 0 - 0 . 8 - 0 . 4 - 0 . 1 + 0.0 
. 3 - 0 . 7 - 2 . 0 - 2 . 3 - 0 . 3 - 0 . 8 - 0 . 1 + 0.1 
. 1 - 0 . 6 - 1 . 9 - 2 . 1 + 0 . 1 - 0 . 5 + 0 . 6 + 0.2 
. 3 - 0 . 3 - 1 . 0 - 0 . 8 - 0 . 2 - 0 . 2 + 0 . 6 + 0.5 
. 3 - 0 . 2 - 0 . 2 - 0 . 3 - 0 . 5 + 1 . 7 + 0 . 8 + 0.5 
. 1 + 0 . 1 + 0 . 5 - 0 . 1 + 0 . 5 + 1 . 4 + 1 . 1 + 0. 3 
. 3 - 0 . 2 - 0 . 1 - 0 . 7 + 0 . 3 + 1 . 1 + 1 . 4 + 0. 4 
. 3 - 0 . 3 - 0 . 4 - 1 . 1 - 0. 0 + 1 . 0 + 1 . 0 + 0. 1 
. 3 - 0 . 5 - 0 . 7 - 0 . 9 - 0 . 2 + 1 . 1 + 0 . 7 - 0.0 
. 5 - 0 . 7 - 1 . 0 - 1 . 1 - 0 . 3 + 0 . 5 - 0 . 2 - 0. 5 
. 7 - 0 . 8 - 1 . 0 - 1 . 4 - 0 . 2 + 0 . 4 - 0 . 9 - 1. 9 
• 8 - 1 . 0 - 1 . 2 - 1 . 2 - 0 . 9 - 1 . 3 - 2 . 0 - 1.8 
. 5 - 0 . 9 - 1 . 3 - 1 . 4 - 1 . 5 - 1 . 7 - 1 . 7 - 1. 9 
. 7 - 1 . 1 - 1 . 5 - 1 . 5 - 1 . 6 - 1 . 9 - 1 . 8 - 1. 9 
. 7 - 0 . 9 - 1 . 2 - 1 . 6 - 1 . 9 - 2 . 1 - 2 . 3 - 2.1 
. 8 - 0 . 8 - 0 . 9 - 2 . 1 - 2 . 9 - 3 . 5 - 0 . 1 - 0. 9 
. 0 - 0 . 7 - 0 . 6 - 1 . 1 - 0 . 7 + 1 . 0 + 0 . 3 - 0.2 
. 6 - 0 . 5 - 0 . 6 - 0 . 4 - 0 . 0 + 0 . 9 + 0 . 8 + 0 . 2 
. 4 - 0 . 3 - 0 . 5 + 0 . 3 + 0 . 9 + 0 . 9 + 0 . 6 + 0. 2 
. 3 - 0 . 1 + 0 . 4 + 0 . 7 + 0 . 6 + 0 . 6 + 0 . 6 + 0. 3 
. 3 + 0 . 1 + 0 . 8 + 0 . 4 + 0 . 5 + 0 . 6 + 0 . 5 + 0 . 4 
. 2 + 0 . 4 + 0 . 5 + 0 . r + 0 . 5 + 0 . 5 + 0 . 5 + 0.5 
. 2 + 0 . 2 + 0 . 7 + 0 . 5 + 0 . 5 + 0 . 6 + 0 . 6 + 0.5 
. 1 + 0 . 7 + 0 . 4 + 0 . 5 + 0 . 5 + 0 . 6 + 0 . 5 + 0 . 5 
. 1 + 0 . 5 + 0 . 4 + 0 . 6 + 0 . 6 + 0 . 6 + 0 . 6 + 0. 5 
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ABSTRACT 


Simulations of the gravity data to be expected from a Lunar Polar 
Orbiter spacecraft utilizing either a Doppler velocity tracking system or 
a gravity gradiometer instrument system are generated using a point 
mass model that gives an excellent representation of the types of 
gravity anomalies to be found on the moon. If the state of the art in 
instrumentation of both systems remain at the level of ±1 mm/ sec at 
10 sec integration time for the Doppler velocity system accuracy and 
at ±1 Eotvos^ at 10 sec integration time for the gravity gradiometer 
system accuracy, inspection of the simulations indicates that a gravity 
gradiometer system will give science data with better resolution and 
higher amplitude -to -measurement noise ratio than the Doppler velocity 
system at altitudes below 100 km. The error model used in the study 
is one where the system errors are assumed to be dominated by the 
point measurement noise and data quantization noise. The effects of 
other, more controllable, systematic error sources are not considered 
in this simplified analysis. For example, both systems will be 
affected by errors in LPO orbital altitude and position knowledge, space- 
craft maneuvers, and data reduction errors. In addition, a Doppler 
tracking system will be sensitive to errors. produced by spacecraft 
acceleration (from outgassing or solar pressure) and poor relative 
position of the LPO, Relay Satellite and ground tracking station, while a 
gravity gradiometer system will be sensitive to errors from spacecraft 
attitude and angular rates. These preliminary study results now need 
to be verified by a more complete error analysis in which all the 
uncertainties of the data gathering process are formally mapped into 
uncertainties in the resulting gravity maps. 

1 Eotvos = 10' 9 sec*^ = lO’ 9 gal/cm = 0.1 mgal/km * 10*^ g's/m. 



I. 


INTRODUCTION 


To date, the information that we have on the gravity field varia- 
tions of the moon (and the rest of the bodies in the solar system) has 
come from measurements of the changes in velocity of a spacecraft 
passing by the body, usually by measurement of the Doppler frequency 
shifts in the spacecraft transmitter carrier frequency. This technique 
requires the constant use of a ground tracking station during data collec- 
tion and results in a measurement of one component of the gravitational 
force field (that component along the line of sight to the tracking station). 
However, since nearly constant tracking was usually needed for other 
purposes and the modifications to the spacecraft transmitter to stabilize 
it for use in a Doppler tracking system were minor, the gravity data was 
essentially "free, " except for the costs of processing the data. 

In the coming years, it is proposed to use a Lunar Polar Orbiter 
spacecraft to carry out a complete, high resolution survey of the moon 
from a low polar orbit. One of the more important physical parameters 
to be obtained from the LPO mission is a high resolution gravity map, 
for gravity is one of the few means for obtaining subsurface geophysical 
information. 

The presently used gravity measurement technique using earth 
tracking stations will only give gravity information on the front side of 
the moon, thus, in order to obtain a complete gravity survey of both 
frontside and backside using the Doppler velocity tracking technique, it will 
be necessary to add an 86 kg Relay Satellite to the mission. In the present 
preliminary mission design, ^ the Relay Satellite will be in a high 
equatorial orbit and the LPO in a near polar orbit. The orbit configura- 
tions are such that only 60% of the LPO backside passes can be tracked 
with the Relay. Thus, to obtain a complete, high resolution (15 to 
30 km), gravity survey of both frontside and backside will require 
months of ground station tracking time carefully chosen ahead of time to 
cover the backside gaps caused by the occasional poor position of the 
Relay Satellite. The resultant gravity data will consist of only vertical 
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and east-west gravity components near the equator, and only horizontal 
components near the pole. 

An alternate technique exists for the measurement of the lunar 
gravity field. This involves measuring the gradients of the gravity with 
instrumentation onboard the spacecraft. There are four instrumental 
approaches to the instrumentation and the state of development of three 

12 3] 

of them ' is such that they can be considered as candidate instruments 
for the LPO. Despite some current misconceptions,^ all the gradiom- 
eter instruments can be placed anywhere on the spacecraft, will operate 
satisfactorily at the spacecraft attitude rates presently planned, and will 
have minor spacecraft interface problems. If the gradiometer system 
used is a complete one and measures all the components of the gravity 
gradient tensor, then accurate knowledge of the spacecraft attitude or 
altitude is not needed since it can be extracted from the gradient data. 

The purpose of this paper is not to debate the relative engineer- 
ing merits of the various instrumentation techniques for measurement 
of the gravity field of the moon, but to address the question: 

Suppose both a Doppler velocity tracking system with an 
accuracy of ±1 mm/sec at 10 sec and a gravity gradiometer 
system with an accuracy of ±1 Eotvos at 10 sec were 
equally available, what would be the scientific return 
we could expect from the two systems and how do they 
compare? 

We will address this question by carrying out a simulation of a Doppler 
velocity system response and a gravity gradiometer system response to 
a point mass model that produces a gravity field that closely resembles 
an actual lunar gravity field. In order to keep this paper within bounds, 
we will use a simplified error model where the total errors are assumed 
to be dominated by the point measurement (instrumentation) noise and 
the data quantization noise. A real system noise level will also have 
contributions from other sources. Both systems will be affected by 
errors in LPO orbital altitude and position knowledge, spacecraft 
maneuvers, and data reduction errors. In addition, a Doppler tracking 



system will be sensitive to errors produced by spacecraft acceleration 
(from outgassing or solar pressure) and poor relative position of the 
LPO, Relay Satellite and ground tracking station, while a gravity 
gradiometer system will be sensitive to errors from spacecraft attitude 
and angular rates. These error sources will typically produce errors 
with very long (1 to 3 per orbit) or very short (< 10 sec) period charac- 
teristics and should be separable from the desired lunar data (5 to 
500/orbit). 

However, in the sense that ^hese error sources have not been 
modeled, this is a preliminary s*udy that needs to be verified by a more 
complete error analysis in which all the uncertainties of the data 
gathering process are formally mapped into uncertainties in the reduced 
gravity data. 
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II. 


THE LUNAR MODEL 


Since one purpose of the study is to investigate the high resolution 
capability of the gravity measurement systems, it was desirable to 
generate a lunar gravity model that had high resolution features in it 
and at the same time was a reasonable model for the kind of gravity 
variations we expect to see above the lunar surface. 

Fortunately, lunar gravity data of high resolution and dense 

coverage (and therefore relatively high accuracy) were obtained over a 

small section of the lunar frontside during the Apollo 16 mission. In 

the paper "Lunar Gravity via the Apollo 15 and 16 Subsatellites, " by 

[4] 

Sjogren, Wimberly and Wollenhaupt there is a section of Apollo 16 
subsatellite data from -50 to +30 degrees longitude and +4 to +9 degrees 
latitude where the Apollo 16 Subsatellite was well below 2 0 km from the 
lunar surface. Since the Apollo 16 subsatellite had an inclination of 
only 10 degrees, the track spacing was very close, and even with track- 
ing limited to every third orbit, the data was obtained at 15 km intervals. 
In all, about 58 orbits and 10,000 20 sec observations were available. 

The system accuracy quoted was 1 mm/sec for each 20 sec observation 
point. The Doppler data was smoothed with a 1 min long spline fit and 
then the spline fit analytically differentiated to obtain an estimate of the 
gravity at approximately 15 km intervals. The resultant contour plot 
is the line-of-sight acceleration of the spacecraft at the spacecraft 
altitude. For the central portion of the data this turns out to be very 
close to the vertical component of the gravity vector at about 15 km 
altitude . 

Because this available data s^t had a number of interesting high 
resolution features it was decided to utilize this data as a starting point 
for the generation of a lunar gravity model. Unf o r t unate ly , the data set 
itself, consisting of only the vertical gravity component at only one 
altitude is a little difficult to utilize if we wish to predict the total gravity 
vector, the Doppler velocity, and gravity gradients in different directions 
than vertical and at other altitudes than 15 km. 



III. 


LUNAR POINT MASS MODEL 



To provide the basic gravity anomaly mouel for our simulations, 
we generated a lunar point mass model consisting of 28 large point 
masses buried at various points under the lunar surface and a grid of 
smaller point masses spaced every 15 km along the lunar surface. 

These masses were adjusted by an iterative process until the vertical 
gravity field at 15 km altitude calculated from the point mass model 
had a reasonable resemblance to the Apollo 16 data in the region from 
-22° to +27° longitude and +2° to +9° latitude. Although an excellent fit 
was obtained (no point more than 2 milligals off and the rms variation 
less than 0.6 mgal), the fit was not important to the simulation effort 
since all we really needed to start with was a lunar gravity anomaly 
model that was reasonable in amplitude and harmonic variation for the 
moon, and was not so unlike the real lunar variations as to cause ques- 
tions about the validity of the simulation results. In the same vein, 
although we have found a mass point distribution that does a reasonable 
job of reproducing the vertical gravity field variations, it is by no means 
the only mass point distribution, and its nonphysical nature makes it cf 
little value for interpretation of the lunar features, which should most 
likely be modeled with rings and discs. 

As a check on the validity of the point mass model as well as the 
cutoff points for the computer algorithm, we calculated the horizontal 
gravity gradients (T and r ) at 15 km altitude using the point mass 
model and compared it with the same gradients calculated by differen- 
tiating the original gravity data set, and again obtained excellent agree- 
ment. Other checks for internal consistency were also made, such as 
calculating the horizontal gradients of vertical gravity at 30 and 60 km 
altitude and checking against the horizontal differences of the vertical 
gravity calculated at that altitude. (This check showed that a 9 x 9 array 
of surface mass points contributing to each gravity point were sufficient 
for internal consistency. ) Once the mass point model had been generated 
and checked for consistency, we could then use it to calculate any 
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component of any quantity (three components of the gravity and doppler 
velocity vectors, and rine components of the gravity gradient tensor) 
at any altitude higher than the original data set (~ 15 km) with a high 
degree of confidence that it would give a reasonable approximation to the 
kinds of gravity anomaly variations and measurement system outputs we 
would expect to find over the moon. 
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IV. 


GRAVITY FIELD OF THE LUNAR POINT MASS MODEL 


As a starting point, we first calculated the vertical component of 
gravity from the mass point model at three different altitudes, 15, 30, 
ar A 60 km. These are shown in Fig. 1. The original lunar data set is 
at the bottom for comparison. These were all contoured at 10 mgal 
intervals. With increasing altitude, we see the expected decrease in 
resolution. This is most evident in the two lobed depression with two 
adjacent subsidiary small depressions in the region around -2° longitude. 
The small depressions are seen at 15 km altitude and disappear at 30 km 
altitude, then at 60 km altitude the two strong depressions have merged 
into a single, slightly oval region. 

Since there is no instrument that can measure gravity directly 
in the free fall environment of an orbiting spacecraft, these curves do 
not represent the output of any particular instrument, but only give us a 
feeling for what we can expect from the reduction of the data obtained 
from systems that measure the integral or derivative of the gravity 
field. 
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V. 


DOPPLER VELOCITY DIFFERENCE FROM THE POINT 
MASS MODEL 

We next calculated the effect of the point mass model on the 
velocity of an orbiting spacecraft. A spacecraft in low polar orbit 
around the moon has a velocity in the near north (south) or x ( -x) 
direction of approximately 1, 7 km/ sec. If we wish to obtain high 
resolution data (1/2°, 30 km wavelength or 15 km resolution), then the 
spacecraft will pass over the 15 km resolution cell in approximately 
8.8 sec. The present state of the art in Doppler velocity tracking is an 
accuracy of ±1 mm/sec for a 10 sec averaging time. Thus, data taken 
every resolution cell will have a comparable accuracy. Until more 
definitive information on the exact performance of an LPO tracking 
system is available, we will use this noise val^e of ±1 mm/sec per 
15 km resolution cell in our calculations. 

Now, the exact velocity of the spacecraft at any point will consist 
primarily of the dominant orbital motion plus small perturbations that 
are the integral of all the gravity anomaly forces that have acted on the 
spacecraft. If we were to calculate the effect of the point mass model on 
the velocity of the spacecraft and plot the resultant velocity directly, we 
would obtain a plot with smoothly varying, gradually sloping contours 
(see Fig. 2(a)). To obtain a measure of the gravity field from this 
Doppler velocity plot, it would be necessary to differentiate the velocity 
to find the velocity change as a function of time (see Fig. 2(b)). With 
real data, this is normally done not by subtracting adjacent velocity 
measurements and dividing by time as we did to get Fig. 2(b) from 
Fig. 2(a), but by fitting a cuoic spline fit curve to a set of points 
(typically 6) about the point of interest, and then analytically differen- 
tiating the spline fit curve to obtain the derivative of the Doppler velocity, 
which is the acceleration. This procedure has the advantage that it 
allows for a smooth derivative to be taken and averages out the effects 
of the noise at a single measurement point, but the disadvantage that it 
also smooths out real gravity variations of high frequency, thus lowering 
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the resolution of the processed data. (See example in Fig. 5 of Sjogren, 
et al.,^^ where a cubic spline fit was found inadequate to represent 
high resolution data. ) We do not use this technique in generating our 
plots, as we wish tc explore the high resolution aspects of both the 
gravity signal and the system noise. 

In the plots in Fig. 2(b) and Fig. 3(a), (b), and (c), what is plotted 
is not the Doppler velocity, but the change in the Doppler velocity that 
would be induced in one resolution cell (15 km) by the gravity in that 
resolution cell. Since the time spent by the spacecraft in passing through 
that resolution cell is just 


dt 


15 km 

1 . 7 km/ sec 


8. 8 sec 


( 1 ) 


then the change in Doppler velocity is just 

dv z = g z dt 

which comes out to be 0.88 mm/sec velocity change for 10 mgal of 
gravity acting for 8.8 sec. Thus these plots of vertical Doppler velocity 
difference for a constant resolution time of 8.8 sec are directly 
relatable to the vertical gravity field. The plots of the vertical Doppler 
velocity change are plotted in Fig. 3 for 15, 30, and 60 km 
altitude. At the bottom is the vertical gravity at 15 km for comparison. 
The 15 and 30 km plots are contoured at 5 mm/sec change per resolution 
time or 5 times the single point noise value and the 60 km plot is con- 
toured at 2 mm/sec change per resolution time or twice the single 
point noise value. The 15 km plot could have been contoured at 
0. 88 mm/sec change per resolution time and would have looked identical 
to the vertical gravity plot. It was contoured at 5 mm/sec change per 
resolution time so that we may later compare it with the gravity 
gradiometer plots and other plots with contour values of comparable 
amplitude -to-measurement-noise ratios. One feature to notice on the 
Doppler velocity change plots is that the amplitude of the larger, low 
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frequency, low resolution peaks such as the Lamont peak at +23 longitude 
and the peak at -18° longitude fall off relatively slowly with altitude. A 
factor of four in altitude between 15 km and 60 km only lowers the peak 
by a factor of 2. 5. 
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VI. 


GRAVITY GRADIENT OF THE POINT MASS MODEL 


To examine the performance of a gravity gradient measurement 
system on an LPO, we next calculated various components and combina- 
tions of components of the gradient of the gravity field. This was done 
by calculating at each measurement point the contribution to the total 
gravity gradient field from each of the 28 large mass points and the 
array of 81 smaller mass points in a 9 x 9 square array around the 
mass point below the measurement point. (This cutoff at ±4 points on 
either side of the measurement point was found to give negligible error 
at both 15 and 60 km altitude. ) Since the gravity gradient field is a 
tensor field, the equations for rotating the tensor field from the point 
mass coordinate frame to the calculation point measurement frame 
involves the double application of the rotation matrix. Once each con- 
tribution from each mass point has been rotated to the same calculation 
point reference frame, the amplitudes can be added. (For those 
interested, the derivation is in the appendix.) 

The first set of gravity gradient calculations was the calculation 
of the vertical gradient of the vertical gravity 


zz 


dz 


which for a point mass M at distance R is 


2GM 


zz 


(3) 


(4) 


Because this is a component of a tensor instead of a vector, the change 
of the amplitude of this component of the gravity gradient tensor with 
horizontal separation from the mass point 



is different than the variation in the amplitude of the vertical component 
of the gravity vector with change in horizontal separation 


g z (x) = 


GM 

(z +x ) 


2 2 ^ 
(a +x ) 


( 6 ) 


This is shown graphically in Fig. 1. Note that the vertical gradient of 

vertical gravity (r ) peaks to a maximum over the mass point, just as 

zz 

the vertical gravity (g ) does (see Fig. 4), but drops off more rapidly 
, z 

with horizontal separation distance and becomes slightly negative at 
large separations. Despite this difference, we see there is a great deal 
of similarity in the general shape of the two curves, and one would 
expect that a contour map of the vertical gradient of vertical gravity 
would look very similar to a contour map of the vertical gravity. How- 
ever, we might expect to see higher resolution data in the gradient map 
because the more rapid fall-off of the point function of the gradient. 

In Fig. 5 we present contour maps of the vertical gradient of the 
vertical gravity calculated from the point mass model for 15, 30, and 
60 km altitude. The vertical gravity contour map at 15 km is placed at 

the bottom for comparison. The units of the gradient plots are in 
* 

Eotvos. The present state of the art in gravity gradient instrumentation 
is a noise level of ±1 Eotvos for a 10 sec integration time. (All three 
instrument concepts under development hav£ demonstrated performance 
at close to these levels in their laboratory physical feasibility models 
and all three engineering models presently being assembled have this 
noise level as their goal, ) Since the time for a spacecraft to pass over 
a 15 km resolution cell is 8. 8 sec, the data point in each resolution 
cell can be expefcted to have an accuracy comparable to this. Until more 


*(1 Eotvos = 10^ gal/cm = 10”^ sec” 2 = 0. 1 mgal/km * 10”^ g's/m.) 
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Fig. 4. Strength of vertical component of field 
as a function of horizontal separation. 
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definitive information of the exact performan e of an LPO gravity 
gradiometer system is available, we will assume a system noise level 
of ±1 Eotvos per 15 km resolution cell in our calculations. 

The contours on the gravity gradient maps are 5 Eotvos for the 
15 and 30 km plots (five times the single point noise level) and 2 Eotvos 
for the 60 km contours, (twice the single point noise level), thus the 
contour interval -to-single point noise level (a rough S/N ratio) is the 
same for these contour maps as for the Doppler velocity change con- 
tour maps in Fig. 2 

There are two features to look for; one is that the higher resolu- 
tion features are more readily resolved with the gravity gradient plots, 
this is especially evident in the features near -2° longitude. The other 
is that the larger, lower frequency, lower resolution features fall off 
more rapidly with altitude. For the Lamont peak at +23° longitude and 
the peak at -18° longitude, the factor of 4 in altitude from 15 km to 60 km 
causes a decrease in amplitude of a factor of 4, whereas in the Doppler 
velocity plots the decrease was only a factor of 2. 5. Nevertheless, 
despite the more rapid decrease, the signal -to -noise of those peaks at 
60 km altitude in the gravity gradient plots (8 Eotvos) is twice as good 
as the signal to noise of the same peaks in the Doppler velocity differ- 
ence plots (4 mm/sec change per resolution time). 

One of the more startling features of the gravity gradient plots 
is that the plot of the vertical gradient of vertical gravity at 15 km 
altitude almost seems to show more structure than is evident to the 
eye when looking at the original set of vertical gravity data. This may 
raise some questions as to the validity or reality of the high resolution 
features in the 15 km gravity gradient data. 

First, we are fairly sure that the simple computational algorithms 
used are correct because of our previous tests of internal self-consistency. 
Second, although our contouring routine leaves something to be desired, 
the contours presented were rounded off by eye and therefore are 
smoother (have less detail) than the data set would have them. Third, 
although we used a real set of gravity data as a starting point, our whole 
simulation and the comparisons between the two gravity measurements 
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systems are based on their performance in measuring the gravity field 
of the point mass model. Since the gravity gradient plots show more 
detailed structure of the point mass model than the Doppler velocity 
plots, we would expect to see the same advantage when they are mea- 
suring the distributed mass variations of the moon. 

Fourth, we generated a second point mass model consisting of the 
original 28 large masses, plus an array of smaller mass points, again 
spaced at 15 km intervals, but now placed 15 km below the lunar surface 
rather than on the surface. The convergence was not as good this time, 
and the milligal variations between the point mass model vertical 
gravity and the original lunar vertical gravity ranged as high as 8 mgals 
with an rms variations of 1, 1 mgals. The plot of vertical gradient of 
vertical gravity from this point mass model at 15 km altitude (30 km 
from the mass points) has nearly all the high resolution features of the 
plot from the surface point mass model, and only differs significantly 
at the large peaks (where the convergence was poorest). This indicates 
that the detail seen in the gradient plots is_ in the gravity plots (although 
difficult to see with the eye) and is not an artifact caused by the close- 
ness of the point mass model. 

Now, the vertical gradient of the vertical gravity was chosen 
for the previous plots because the contour plots of this quantity are 
similar in shape to the vertical gravity contour plots and they can be 
easily compared by eye. Of the various possible instrument mechaniza- 
tions for the gravity gradient experiment, most do not give the vertical 
gradient of vertical gravity directly, but instead give combinations of the 
components of the gravity gradient tensor. If two or three orthogonal 
instruments are used (two for certain types of instruments under 
certain conditions) then all the gravity gradient components (including 

T ) can be calculated from the instrument outputs. However, as we 
z z 

shall show, even one instrument, with its slightly ambiguous combina- 
tion output, can give highly useful information. One type of gravity 
gradiometer instrument produces two outputs 



if oriented in the vertical -east plane or 


l 


t 


• r zz - r xx> 


and 


zx 


( 8 ) 


if oriented in the vertical -north plane (the orbital plane of a polar orbit). 
Now, the vertical gradient of the vertical gravity (r ) and the horizontal 
gradient of the horizontal gravity (T xx or r^y) are closely related to each 
other for most types of mass distributions. For a point mass, the 
relationship is simple and direct 

r = -2r = -2r = 

zz xx yy 3 

whereas for more complex features the relationships are only approxi- 
mate. (The relationship is worst for linear features where the hori- 
zontal gradient along the feature is zero. ) Thus a contour plot of just the 
instrument output 


V r zz 


■ r xx> 


4 


ZZ 


(10) 


can be expected to have nearly all of the features of a contour plot of 
vertical gradient of vertical gravity with some slight distortion. This 
can be seen in the plots in Fig. 6. These are all at 30 km altitude and 
have been contoured at 2 Eotvos (twice single -point-noise -level) to bring 
out the detail. As can be seen by inspection and comparison with the 
vertical gravity plot at 30 km, a simple contour plot of the instrument 
output can provide high resolution information that can be used for 
geophysical interpretation. 

Now the plots of gradiometer instrument output are in units and 
parameters and mixtures of components that geophysicists are not used 
to working with. There are two approaches to alleviating the problem. 
One is to develop algorithms that will remove the ellipsoidal distortion 
to obtain the vertical gradient and others to obtain the vertical gravity. 
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However, with a little retraining, geophysicists may find that these 
plots are just as easy to work with as vertical gravity plots and more 
valuable because of their high resolution content. 

At the present time, when a geophysicist looks at a contour 
gravity map, he compares the features with a topographical map, and 
then using experience and reasonable assumptions, constructs a mass 
model of layers, scarps, discs, rings and other geometrical features. 
He then uses a computer to calculate the gravity from these features 
and compares it with the real data set, and readjusts the model until 
he gets a good fit. The procedure with gravity gradient data could be 
the same, except that the computer program is changed to calculate 
th^ gradient components (or combinations of components) instead of the 
gravity component. The resultant mass distribution should be the 
same for eitner procedure. 
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VII. GRAVITY PLOTS FROM GRADIENT PLOTS 1 


t 


\ 


[ 6 ] 


If the gravity anomalies are desired in the form of gravity 
contour plots, then it is possible to obtain these in a simple manner 
from the gravity gradients by integration of the gradient components. 

The time variation of the components of the lunar gravity 
vector (g) arc given by the product of gravity gradient tensor (r) and the 
vehicle velocity vector (v) 


g 


r . v 


(ID 


Equation (11) can be written component by component in a north, west, 
up (x,y,z) coordinate frame; 


ag x ag x ag x 

g = ~ — v + — — v + ~ — v (north component) 
®x dx x 9y y c)z z r 



ag 

__I V + 

ay y 



(west component) 


ag z ag ag z 

K = a^" v x + aT v y + aT v z (verticalcom P° nent) {12) 


where v^, v^, and are components of the velocity, v. Because the 
gravity gradient tensor is symmetrical and the Laplacian of the gravity 
potential is zero the following relations apply; 


9g x 98 v 

mm rnimmmmmJLmm 

9y ~ 9x 


ag„ ag z ag v ag z 

— — = — — JL s — — 113) 

dz dx ’ dz dy 
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dz 


o 


( 14 ) 



d lx 

dy 


Note that the gravity vector can be computed from eqs. (12) by 
integration with time if the gravity gradients and vehicle velocity are 
known. Note also that if vehicle velocity is constrained so that certain 
components are zero (e.g. , v^ = 0, v^ = 0), then it is unnecessary to 
measure all gravity gradients, because some terms drop out of these 
equations. 

For the case of the LPO spacecraft in a polar orbit, the major 
component of the spacecraft velocity vector is in the north (x) direction, 
while the west (y) and up (z) components of the velocity vector are 
much smaller. For illustration, let us assume that: 


v = 1.7 km /sec 

x 


v = 0 

y 


v = o 
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(15) 


then the eq. (12) reduce to: 
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Thus by utilizing the fact that the velocity of the* vehicle is uniform, we 
can see that from cqs. (16) that it is only necessary to measure three of 
the five components of the gravity gradient tensor to obtain all three 
components of the gravity field vector. If the spacecraft has two 
rotating gravity gradiomete rs , one oriented with its spin axis vertical 
to measure the gradient components in the (x, y) plane tangent to the 
surface of the moon, and the 1 other with its spin axis along the orbit 
axis to measure the gravity gradients in the (x, z) plane of the orbit, 
then the sine and cosine outputs of the sensors will contain the following 
gravity gradient components: 


xy 


zx 


3g x 


8x 

dy 


dg* 

dz 

dx 

of the 

sensors 


xy 


zx 


= 2 


= 2 


de 

ax 

3g z 

9x 


(17) 


gravity gradient components directly. The remaining gravity gradient 
term can be obtained by combining the two sine outputs and using the 
fact that the Laplacian of the gravity potential must be zero (eq. 14) 



-S 


zx 


ffx l^x . 

3x ” V 3x dy 




(18) 


Thus, by utilizing the tact that the LPO velocity vector is nearly entirely 
in the lunar north direction, we can obtain a good estimate of all the 
components of the lunar gravity field vector with the output data from 
only two rotating gravity g radiometers: 
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(19) 


Alternatively, we can obtain a vertical gravity contour map from 
just one of the outputs (r ) of one gradiometer by integration: 


‘/"sr ' "k 61 


■/ 


dx 


dx 


( 20 ) 


If the orbit were truly polar, then the integration would start as the LPO 
passed over the pole. A nominal value for the gravity at the pole would 
be assumed as the constant of integration, and then the gradient varia- 
tions from that point would be used to calculate the gravity variations 
from the nominal value. In practice, the orbit will not be truly polar 
and a fit would have to be made to the circle of closest approach to the 
pole to obtain starting values for each orbit. 

In Fig. 7(a) we have presented a contour plot of the north varia- 
tion of the vertical gravity (r ) as calculated from the point mass 

zx 

model at an altitude of 15 km. Since V goes to zero over a point 
mass, the contours seem to have no relation to the vertical gravity. 
However, if we assume that integration has been taking place in the 
data along the orbit from the pole, then we will know the gravity field 
along the top edge of the region that we are plotting. If we take this 
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single line of vertical gravity data and then add or subtract an amount 
given by 


= r 2s d x BD 

(where dx is 15 km), then for every change of 1 Eotvos in 1 5 km we get 
a change in the vertical gravity of 1.5 mgals. This integrated set of data 
was contoured and plotted below the gradient plot (Fig. 7(b)) and is seen 
to be very close to the vertical gravity at 15 km calculated directly 
from the point mass model (Fig. 7(c)). Since the noise on the gradient 
data is ±1 Eotvos, the noise of the gravity contours is ±1.5 mgals. 



10 mgal 10 mgal 5 EOTVOS 

CONTOURS CONTOURS CONTOURS 



15 km ALTITUDE 


4662-14 




VIII. EFFECTS OF SYSTEM NOISE ON CONTOURS 


It is difficult to give a direct and simple relationship that will 
describe how the noise of the measurement system affects the contours 
used in a plot. The measurement noise is a random variation of a 
single point measurement. In our simulations here, we assume a mea- 
surement every 15 km in both the x and y directions. The LPO space- 
craft will be in a near polar orbit with a velocity in the north-to -south 
or south-to -north direction of about 1. 7 km/sec, thus we are assuming 
a measurement every 8.8 sec. Naturally the measurements along the 
track, being sequential in time and at nearly the same altitude, will have 
more correlation than those across the track, which will come on 
different orbits. 

To get a feeling for the effects of measurement noise on the 
contours, we ignored the correlation effects and assumed that each mea- 
surement point was an independent measurement. We then added 
gausslan random noise onto each calculated point for both the Doppler 
and gradient at 30 km altitude. The resultant noisy data caused havoc 
with our contour plotting routine, so it was necessary to apply some 
adjacent point smoothie to *ake out isolated noise peaks. This smooth- 
ing effectively reduced our rms noise level from ±1 Eotvos and ±1 mm/ sec 
to ±0.6 Eotvos and ±0.6 mm/sec. (This was of course accompanied by 
a degrading of the resolution from 15 km to about 30 km, but there was 
probably not much 15 km resolution data left at 30 km altitude anyway.) 
These smoothed data were then contoured and plotted. This was repeated 
four times with four different sets of additive random noise and the 
contours plotted on top of each other to give an estimate of the ’’noise 
width” of the contours (see Fig. 8). From inspection, it is seen that 
most of the Doppler velocity contours have ’’noise widths” greater than 
15 km, in some areas reaching 45 km, while most of the gradient 
contours have ’’noise widths” significantly less than 15 km. It should 
be emphasized that if we make only one LPO gravity survey at 30 km 
altitude, we will get only one ’’best” contour from that one set of noisy 
data, and that ’ best” contour would lie anywhere within the ’’noise width” 
of the contours described by the four superimposed contours of the plots. 
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IX. 


GRAVITY MEASUREMENT AT HIGHER ALTITUDES 


Most of the simulations presented in this pcper have been at 
relatively low altitudes (15, 30, and 60 km). Although low, these 
altitudes were attained over many regions of the moon with the Apollo 
spacecraft and subsatellites and should be easily attainable with any rea- 
sonable spacecraft and mission design for an LPO mission, especially if 
the mission is planned to allow for initial sum ys at higher altitudes 
(60 km) to obtain a good representation of the gravity field for orbital 
calculation before the spacecraft is moved into a lower orbit for the 
higher resolution surveys. However, for many reasons, the mission 
planners may want the option of operating at higher altitudes. What are 
the comparative merits of the two gravity measurement concepts at these 
altitudes? In Fig. 9 are plots of the vertical Doppler velocity change, 
the vertical gradient of the vertical gravity and the vertical gravity at 
100 km. The contours for the Doppler velocity change are at 1 mm/ sec 
change per resolution time (contour inte rval -to -single -point-noise level 
of 1) and the contours for the gradient are 1 Eotvos (similar signal-to- 
noise). We have plotted the contours directly without noise and feel that 
a direct comparison is valid since any measurement system integration, 
smoothing, spline fits, or other tricks to take advantage of the correla- 
tion between adjacent measurement points to reduce the noise level on 
the contours could be applied equally well to both Doppler velocity and 
gradient data. It seems obvious from inspection that even at ^ km 
altitude the gravity gradient data is superior in both resolution id 
amplitude -to -single -point-noise level. This does not hold true for 
higher altitudes, however. At ZOO km the amplitude variations of the 
Doppler difference plot and the gradient plot have comparable S/N ratios 
(about 1:1) and at altitudes higher than ZOO km, a 1 mm/sec at 10 sec 
Doppler system will outperform a 1 Eotvos at 10 sec gravity gradiometer 
system in exploring the larger anomalies. 
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X. 


SUMMARY 


In summary, a set of simulations of the gravity data to be 
expected from an LPO spacecraft utilizing the state of the art in both 
Doppler velocity tracking systems and gravity gradiometer instrument 
systems has been simulated using a point mass model that gives an 
excellent representation of the types of gravity anomalies to be found 
on the moon. If the state-of-the-art in instrumentation of both systems 
remains at the levels of ±1 mm/sec for 10 sec integration time for the 
Doppler velocity tracking system and ±1 Eotvos for 10 sec integration 
time for the gravity gradiometer system, it is seen by inspection of the 
simulation plots that the gravity gradiometer system will give higher 
resolution data with higher signal-to-measurement noise ratio than the 
Doppler tracking system at all altitudes below 100 km. The simulations 
concentrated on instrumentation errors and did not cover many of the 
systematic error sources involving spacecraft altitude, attitude, angular 
rate, position or accelerations. It is recommended that the simulations 
be redone using a more complete error model and a more rigorous 
error analysis to quantify the differences in performance of the two 
systems. 
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ROTATING COMPONENTS OF THE TENSOR GRAVITY GRADIENT 
FIELD OF A POINT MASS TO A MEASUREMENT POINT COORDINATE 
SYSTEM FROM THE POINT MASS COORDINATE SYSTEM. 

As shown in Fig. A-I, a point mas s at the position (x^ , , 0) 

produces a gravity gradient tensor field at the measurement point 
(0, 0, h ) that is tilted with respect to the x, y, z lunar coordinate sys- 
tem. To determine the gravity gradient tensor components in the 
lunar vertical coordinate system we must determine the rotation 
matrix coupling the two coordinate systems. 

To rotate from lunar coordinates to the coordinate system 
defined by the mass point and the measurement point we first 
rotate the lunar (x, y, z) coordinate system about z until the y axis 
goes through (xj,yj,0). This is a negative 4> rotation that moves y to the 
H axis and x to the £ axis, with the 4 axis remaining identical to the 
z axis. We then rotate the (£ f n f £) coordinate system an angle 0 about 
the £ axis through the point (0, 0, h) to tilt the C, (vertical) axis until it 
intersects (xj,yj,0). The rotation moves the ^ axis to the y 1 axis 
the axis to the z 1 axis while the % axis becomes the x 1 axis. The 
principal axes of the gravity gradient tensor of the mass point are now 
seen to be oriented along the (x'.y'jZ 1 ) coordinate system. To rotate 
the tilted tensor into lunar coordinates, reverse the above. Using the 
general rotation matrix from Goldstein page 109, ' with ip = 0 
(sin ^ = 0, cos 4*= 1). The rotation matrix is: 

cos 4> sin $ 

-cos 9 sin 4> cos 9 cos 4 

sin 0 sin $ -sin 6 cos $ 




H. Goldstein, Classical Mechanics, Addison-Wesley, Cambridge, Mass. 

( 1950 ). 
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and it f s inverse is: 


s:, 1 

j i 



•cos 0 sin 4> 
cos 0 cos 4> 
sin 0 



The gravity gradient tensor at a point (0, 0, h) from a mass M at 
(x, y, 0) is given in the x', y', z' system as 


-1 




GM 
„ 3 


0 0 


-1 


We now want to rotate this to the lunar x,y, z system, so we reverse 
the rotation 


C 0 = sT. r!. s.„ 

k£ ki ij j£ 


to obtain the gravity gradient components in the (x,y,z) lunar 
coordinate system: 


1 

[-1 + 3 sin 2 8 sin 2 ♦ J 
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From the geometry of Fig. A-l we see that: 
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cos 8 



sin 4* 
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sin 0 = 


These can be used in place of the angles in computing the tensor 
components. (Take care to reduce the equations to eliminate p from 
the denominator since p = 0 for a mass point directly under the 
measurement point.) 
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Fig. A-l. 

Rotation of gradient 
tensor. 
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1.0 SCOPE 


This document covers the functional and physical description 
and the functional performance parameters of the Gravity 
Environment Measurement Subsystem (GEMS). In addition, the 

sensor/spacecraft interface characteristics and constraints are 
defined. This document is intended to be a working document with 
the addition of data following selection of the specific 
spacecraft . 

2.0 Applicable Documents 

(1) MJU Spacecraft Description Document 

(2) JPL Document 660-6 

(3) "Rotating Gravity Gradiometer Study", Final Report on 
JPL Contract 954309 under NASA Contract NAS7-100 for 
the period from August 1975 through June 1976, Hughes 
Reseaich Laboratories, Malibu, CA 90265. (The final 
report of which this document is an Appendix.) 


3.0 Functional Description 

3.1 Scientific Objectives 

The objective of the Gravity Experiment is to obtain a high 
resolution, free-air gravity map of the moon. The gravity data 
will be combined with orbital, Doppler tracking, altimetry and 
photographic data to obtain a high resolution global gravity map 
for : 

a. Orbital dynamics predictions 


b. Geodesy 

c. Geophysics 


d. Cartography 

Gravity gradient data may also be combined with )ther data 
to improve post-mission knowledge of spacecraft orbit and 
attitude. 


3.2 Experimental Approach 

The Gravity Environment Measurement System (GEMS) , which 
consists of one or more Rotating Gravity Gradiometer (RGG) 
sensors and supporting electronics, will be used to obtain the 
in-situ lunar gravity field data. A general block diagram of a 
typical GEMS for space application is illustrated in Figure 1. 
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Figure 1. Gravity Environment Measurement System (GEMS) 
for Space Applications 


The system illustrated comprises two Rotating Gravity 
Gradiometers (RGG's) with internal electronics, an external 
electronics unit for eacn RGG, and common power supplies and 
interface units. The LPO spacecraft (S/C) support equipment is 
also shown to indicate the interface requirement. The interface 
is functionally the same regardless of the number of RGG's in the 
system. However, the number of data and command regi&ters would 
vary. 


Each RGG instrument or "sensor" has two output signals. One 
of these is a direct measure of the cross gradient tensor 
component in the plane of the sensor rotation and the other 
signal is a measure of the difference of the two principal 
gradient tensor components in the plane of rotation. The 
"principal gradient difference" outputs of three mutually 
orthogonal sensors are required to determine the principal 
gradient tensor components separately. However, for the LPO 
application, studies indicate that it is possible to obtain 
adequate data for determining subsurface mascon densities with 
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two, or even one RGG sensor. This is discussed in detail in Task 
3 and 5 of document (3) referenced in Section 2 above. 

For simplicity, the following discussions will be based on a 
one-RGG system. 

3.3 Major Functional Elements 

3.3.1 RGG Sensor 


Figure 2 is a photograph of a prototype RGG. The connector 
panels at top left and bottom left and the base against which the 
rule is leaning are temporary items used for laboratory testing. 
In operation, the RGG is mounted by the center flange which is 
about 23 cm in diameter. The RGG is 30 cm long as pictured, but 
will be slightly longer when cooling fins are added at each end. 

3. 3. 1.1 Arm/Pivot Structure 


The basic elements of the RGG sensor are shown diagra- 
matically in Figure 3. Two arms with masses, M, at the ends are 
oriented at 90 degrees with respect to each other. They are 
mechanically coupled by a torsional pivot spring and supported 
within a housing by additional torsional pivot springs. 

3. 3. 1.2 Transducer 


Figure 4 shows a piezoelectric transducer attached to the 
arms. The supports cause the transducer to bend in proportion to 
the differential motion of the arms about the pivot axis. This 
bending action produces an output voltage proportional to the 
gravity gradient field. 
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3. 3. 1.3 Spin Bearings 


The arm/pivot/transducer/inner-case assembly called the 
rotor, which also includes other elements such as the internal 
electronics discussed below, is supported on a pair of precision 
hydrodynamic oil bearings, as diagrammed in Figure 5. 






3. 3. 1.4 Spin Motors 

The rotor is spun within an evacuated outer case called the 
stator by a pair of asynchronous eddy-current drag-cup motors, 
one at each end, as diagrammed in Figure 6. 


S27M 


DRAG CUP 




30 cm 



Figure 6. Addition of Motors, Transformers and Stator 


B-12 




3. 3. 1.5 Speed Sensing 


An eight-slot disc at the base of one drag cup is pictured 
in Figure 7. Optical sensing of the slots, one of which is wider 


5278-19 



Figure 7. RGG Drag Cups, One with Speed Control Slotted Disc 


than the other, plus digital processing and a closed loop servo 
provides the rotor-to-stator position, speed data, and speed 
control. This is covered in more detail under electronics. 
Section 3.3. 


3 3.1.6 RF Transformers 


The sensor has two rotary transformers, also shown in 
Figure 6. One is used to transmit PF power and logic signals to 
the rotor to supply power for the internal electronics and 
control the sensor output mode condition. The other transformer 
is used to transmit gravity gradient and test signals from the 
rotor to external electronics. 
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3.3.2 Electronics 


The GEMS electronics and related elements for space appli- 
cations, as shown in Figure 8, contains three parts. The portion 
which is contained within the RGG is shown at the left within the 
symbolic RGG outline; related electronics, which are external to 
the RGG, are shown in the center; and the spacecraft support 
equipment is at the right. However, for easier comprehension, 
the electronics will be divided into funr .ional portions in the 
following discussions. 

3.3. 2.1 Speed Control Elements 

Figure 9 shows the GEMS functional elements which constitute 
the speed control . A light emitting diode (LED) and photo diode 
sensor (under the box protruding from the sensor case at the 
upper right in Figure 2) detect the eight slots in the slotted 
disc (See Figure 7). This signal is amplified and shaped by the 
photo pickoff providing eight digital pulses per revolution of 
the rotor (3 w pulses) to the digital interface. Here the time 
between octant pulses is determined by counting the output of a 
high precision 10 MHz clock. 

The count for each octant is corrected for irr-gular slot 
spacing and compared with a commanded number representing the 
required speed. The difference between the corrected octant 
count and the command speed number is added to a long-term error 
count forming the speed error value. This digital number is 
converted to an analog signal by the D/A converter and sent to 
the motor driver and 2-channel amplifier to drive the RGG motors 
to a higher or lower speed. This correction is accomplished 8 
times per revolution and maintains the rotor speed to within one 
part in ten mil) ion. RGG speed control timing and flow c K rt are 
shown in Figures 10 and 11. 




Figure 8. GEMS KGG and Electronics for Space Applications 





































Figure 9. GEMS Speed Control Elements 



























3. 3. 2. 2 Power Supply and Therma l 
Control Electronics 

GEMS receives raw 28 VDC power from the S/C. This raw power 
is converted to the required power forms in the external 
electronics or into RF for transmission through one of the rotary 
transformers to the RGG rotor. The RF is rectified and filtered 
to provide DC for the RGG internal electronics. It is assumed 
that raw power switching will be performed by the S/C power 
subsystem. 

Temperature control of the RGG stator begins with a 
commanded temperature operating point received via the digital 
interface, as shown in Figure 12, and set into the temperature 
control unit. Two temperature sensors are used; one measures the 
ambient temperature and the other the RGG case temperature. The 
latter sensor is a special wire wound in a small helical slot on 
the periphery of each half of the stator and a heater wire is 
likewise wound in a larger helical slot in each half interleaved 
with the sensor wire slot (See Figure 2) . The sensors and 
heaters together with the temperature control and power unit 
provide a closed loop servo to maintain the RGG case within 
+ 0.06 C of the commanded operating temperature point. 

3.3.2. 3 Data Output Electronics 

Gravity gradients stimulate the rotating arm/mass/pivot/ 
transducer structure to oscillate with an alternating angular 
motion at the 2 w (approximately 35 Hz) resonant frequency, 
producing an amplitude modulated signal from the piezoelectric 
transducer. Operating at resonance results in (1) a mechanical 
amplification of the gradient induced differential arm motion of 
about 300 (Q approximately equal to 300) and (2) mechanical 
filtering of noise signals before the transducer converts the 
mechanical motion to an electrical signal. This significantly 
increases the signal-to-noise ratio over that which could be 
attained with electronic amplification alone. 

The electrical signal is amplified, first by a balanced pair 
of low-noise field-effect transistors (FET's) in the 
preamplifier, and then by an amplifier with controllable 
normal-gain and low-gain modes (See Figure 13) . Additional 
electronic tuned filtering is also provided. 

The amplitude-modulated signal is converted to a 
frequency-modulated (FM) signal by a 200 KHz mid-frequency 
voltage controlled oscillator (VCO) type modulator. This FM 
signal is then transmitted to the external electronics through a 
rotary transformer. Another rotary RF transformer at the other 
end of the sensor is used for transmitting powefc into the RGG 
rotor. 

























) 


The FM sensor signal from the RGG rotary transformer goes to 
an FM receiver, then to a zero-crossing counter and latch with a 
digital output which is directly related to the voltage from the 
piezoelectric transducer and thus directly related to the sensed 
gravity gradient. The counts for each octant of rotation 
(quadrant of space) are averaged over a selectable period, 
typically 64 rotations or approximately 3.6 seconds. The output 
is fed to the digital interface for output to the spacecraft data 
system. 


3. 3. 2.4 Test and Control Electronics 


GEMS test and control functions related to the RGG are 
illustrated in Figure 14. 

(1) Logic Control ? The RF power control transmitted into 
the rotor is interrupted in a time sequence digitally coded 
signal by the data injector unit. This signal is received and 
decoded by the interrupt detector for use in the logic control of 
the functions discussed in (2) through (6) below. 

(2) Amplifier Gain : This logic control selection provides 
for normal gain or low gain of the sensor output amplifier. Low 
gain is used during sensor check out and test to avoid 
overdriving the amplifier output stages when the bias or signal 
level is higher than normal. 

(3) Voltage Test : An accurate bridge measurement of the RGG 
internal raw or regulated dc voltage can be selected for output 
in place of the gravity gradient signal. Simultaneously, an 
input reference 2 w square wave voltage is switched to the 
voltage test bridge. 

(4) Torque Test : A test and calibration torque can be 
applied to the RGG arms by routing of the input reference 2 w 
square wave to capacitor plates located near the end masses. The 
2 w signal is correctly phased by the digital interface unit 
using the 8 w pulses from the photo pickoff unit as a reference. 
Because the output signal from the transducer is a sum of the 
gravity gradient signal and the calibration output signal 
produced by capacitor torquers, the same preamplifier, amplifier 
and output switch logic selection are used. 

(5) Temperature Test ? Internal RGG temperature measure- 
ments are made by logic switching of the input 2 w square wave to 
the temperature test bridge. This along with the input from the 
temperature detector produce an output which is simultaneously 
selected for transmission to the external electronics in place of 
the gravity gradient and torque signal. 
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Figure 14. GEMS Test and Control Electronics 












(6) GEMS Input/Output : Command data to GEMS from the S/C 
computer, or from earth, and science and engineering data from 
GEMS for storage in the S/C recorder or for transmission to earth 
are handled through the digital interface unit's input/output 
(I/O) addressable registers. Clock and sync signals for I/O 
timing reference are received from the S/C through the analog and 
discrete interface unit. I/O is discussed in greater detail in 
Sections 6.4 and 6.5. 

3.4 Operating Modes 

The GEMS operating modes are controlled by command and vary 
only in the characteristics of the data being acquired (science 
or engineering) , the spin speed, sensor operating temperature, 
etc. Therefore, the GEMS has only three power modes: Standby 
(temperature stabilization); Operate; and Off. 

3.5 Preliminary Space GEMS Electronics Subsystem Design 

A first look at the electronics subsystem design for the LPO 
and other space applications was done with the primary purpose of 
determining interface requirements and approximate size and 
weight of the external electronics. It was based on the use of 
standard Hughes space-qualified and proven design and fabrication 
techniques and hardware. In the process of accomplishing this 
task, the functional block diagrams of Figures 15 and 16 were 
generated along with other information included in other parts of 
this report. For specific information on size, weight, power, 
etc, see Section 6.0. 

Figure 15 shows the functional relationships of a number of 
electronic elements which require digital registers and the 
register sizes. Figure 16 shows GEMS-Spacecraf t command and data 
interface requirement including timing and control lines. 

Lists of the command, science and engineering data items 
with the requirements for registers (words) , register sizes 
(bits/word) and data rates are included in Sections 6.4 and 6.5. 




Figure 15. Preliminary Register Functional Block Diagram for Space GEMS 


























4.U Ph¥SICAL LESCblPTIOW 


In the Discussion of GEhS above, the simplif ied 
representation of the EGG was used for ease of functional 
Description. There are seme differences between the simplified 
version ana the actual hardware which arc describee in the 
following paragraphs. 

4.1 Sensor Arms 


bather than the single "bar" type arm depicted in the 
functional Diagrams, each sensor arm is actually a pair of arm 
plates with the end masses bolted between as shown in the 
photograph of Figure 17. This forms a rigid box structure 
resistant to bending under high-g shock and launch conditions. 
It also reduces errors due to anisoelasticity of the arms. 
Cutouts in the arm plates are made to render the arms isoelastic 
and minimize these errors. The arms are interleaved to get the 
centers of mass closer together and for better position of the 
torsion spring supports. 

4.2 Inner Kotor Housing 

The arms with their end masses are enclosed in a case 
comprising a circular center plate shown in Figure 17 and two end 
bells shown in Figure lb. The bells are bolted to the center 
plate and to the bases at the ends of the pivots (one of which is 
shown at the top of Figure 17) forming a very rugged dual box 
structure. 

The relationship between these parts and others discussed 
below is shown in Figure 19. 
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Figure 17. RGG/End Mass/Center Flate Structure 


4.3 Pivot Springs 

Because of the dual-plate, interleaved arm design and the 
center-plate, dual-bell inter rotor housing design, there are six 
torsional pivot springs as shown in Figure 19. Two of these are 
between arm plates of different arms and tend to keep the arms 
orthogonal but permit them to rotate slightly relative to one 
another. The other four pivots keep the arms fixed to the inner 
rotor case except for slight rotation of the arms relative to the 
case. 
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Figure 18. RGG Arm/Mass/Center-Plate Assembly with Bottom 
Inner Rotor Bell in Place and Top One Held 
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Figure 19. Relationship of RGG Structural Parts 


Five of the six pivots are visible in Figure 20. They vary 
from about 3 to 4 mm in diameter and 2.7 to 2.8 mm in length 
depending on position. 

4.4 Arm/Pivot/Center Plate Fabrication 

All of the structure comprising the four arm plates, center 
plate, six pivots and two end bases are machined from a single 
forged aluminum billet and remain a single Piece of «etal 
throughout fabrication. This avoids a very difficult assembly 
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Figure 20. RGG Arm/Pivot/Center Plate and End Mass Assembly 
Showing Top Five of the Six Pivots 


and fastening process, which would be required if these 
functional elements were machined as separate pieces. Figure 21 
shows a section of a forged aluminum billet at the upper right, a 
rough machined part at the left and a final machined part at the 
lower center. 
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4.5 


Transducers 


Instead of a single transducer as depicted in the functional 
sketches, there are two. Each one passes through a hole in an 
enlarged portion of the two pivot sections connecting arm plates 
of opposite arms. These rectangular holes are visible in 
Figure 20. Thus, one transducer is attached between the two 
upper plates and the other between the two lower plates. 
Relative rotation of the arms produced primarily by gravity 
gradients causes the transducers to bend and produce an output 
signal, while rotation of the arms in unison relative to the case 
does not bend the transducers or produce an output. Thus, 
vibration of the instrument caused by environmental disturbances 
do not produce first order errors in the output. Slight errors 
can be induced by coupling of the vibrations with residual 
imperfections in the instrument. Adjustments are provided to 
reduce these errors to levels below the design sensitivity of the 
sensor . 


4.6 Outer Rotor 

The inner rotor case supporting the arm/end mass/pivot 
structure, defines the torsional pivot axis. To minimize one of 
the error sources reguircs adjustment of the pivot axis for 
alignment with the rotor spin axis. Thus, an outer rotor case 
comprising a pair of end bells similar to the inner rotor end 
belis but slightly larger, are also attached at the center plate 
and support the rotating halves of the two bearings. The top one 
is diagrammed in Figure 19. the two motor drag cups and the 
slotted disc are also attached to the outer rotor. The RGG 
sensor with half of the stator removed is shown in the photograph 
of Figure 22. The upper half of the rotor appears protruding 
from the bottom half of the stator. At the top is the upper drag 
cup and around its base is the rotor half of the ring capacitor 
used to couple the 2 w square wave reference and torque driver 
signal to the rotor. The slotted disc shown in Figure 7 is 
attached to the opposite end of the rotor at the base of the 
second drag cup. 


4.7 Stator 

The outer case of the RGG, called the stator, supports the 
rotor by the spin bearings as shown in Figure 19. It contains 
the stationary parts of the two rotary transformers, the 
capacitance coupler and the two spin motors. It also supports 
the optical sensor for speed sensing and provides an ar-tight 
chamber which is evacuated after assembly. The vacuum reduces 
windage on the rotor and thermal conductivity between the stator 
and the rotor. 
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The completely assembled prototype RGG sensor (No. 1) is 
shown in the photograph of Figure 2. See the first paragraph in 
Section 3.3.1 for further information. Heating wires and 
temperature sensing wires can be seen in the helical grooves on 
each stator half. Fins will be added to the case ends for added 
radiative heat dissipation in space applications. See Section 
5.5 for report on the thermal control study. 

4.8 Internal Electronics Packaging 

The internal RGG electronics is located on two printed 
circuit boards one at each end of the rotor between the inner and 
outer rotor bells. The boards for prototype RGG No. 1 are shown 
in Figures 23 and 24. They are slightly saucer shaped and 
attached to the outside of the inner rotor as shown in Figure 19. 
The internal electronics circuit boards for prototype RGG No. 2 
will be flat and attached to the inside of the outer rotor. This 
change will result in mere uniform temperature of the pivots 
which are the most critical temperature sensitive elements in the 
RGG. 


4.9 External Electronics Packaging 

The prototype RGG external electronics comprises mostly 
standard laboratory electronics equipment and special units built 
by Hughes for laboratory testing. However, a preliminary design 
analysis has been made of the external electronics requirements 
for the LPO mission gravity experiment- and other similar gravity 
survey space missions. This analysis indicates that the external 
electronics can be contained in one package 8x15x137 cm. 
Figure 25 is an outline and mounting drawing of the proposed 
package. Design and fabrication would be in accordance with 
standard Hughes space qualified electronics equipment, which has 
proven highly reliable in many space systems. 
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Figure 23 


RGG Internal Electronics Top Circuit Board 









5.0 SPECIAL REQUIREMENTS 


5.1 Sensor Location on Spacecraft 


The rGG and electronics can be placed at any convenient 
location on the S/C. '..here is no requirement to locate the 
sensor at the S/C center-of-qravity (CG) nor to control the 
spacecraft CG. If the sensor is located in close proximity to 
either movable devices (scan platform) or expendables (fuel) 
these effects can be remeveo from the data by modeling or by 
disregarding the data taken curin', periods when large masses arc- 
moving . 


One factor which may affect location of the wGG is power for 
thermal control. Tc keep this near the lew end of the ranee, it 
is desirable tc place the FCG ir. a location which has a low- 
thermal gradient. This is discussed in detail in Section 5.5. 

a.2 Sensor Alignment 

The sensor spin axis alignment, both control and knowlecqe, 
is not critical. Nominal S/C alignment tolerances cf 1 tc 5 
degrees will suffice since alignment errors only cause a nominal 
initial bias which will be removed from the data along with other 
initial biases including the S/C gravity gradient. Actual low- 
frequency GEMS gravity data plus what is previously knev*n about 
the lunar gravitational field ir the low frequencies can be used 
to later calculate the sensor attitude with respect to the mean 
Lunar vertical to better than a millirodian (See Appendix D) . 


5.3 Sensor Attitude Control Rates 


The GEMS experiment uses a "specific force" gradiometec and 
if the S/C is undergoing rotation, then the induced rotational 
velocity field produces a rotational gradient that is sensed by 
the grediometer. The gradient produced by a rotational rate has 
an amplitude equal to the square of the rotational rate ar.d is 
cylindrically uniform about the axis of rotation. The following 
two conditions have been considered. 

(a) If the S/C is gravity gradient stabilized it is rotating 
about the orbit axis at a rate equal to the instantaneous orbital 
angular rate. There will then be a positive rotational gradient 
that is uniform in the plane of the orbit and whose amplitude 
(for a circular orbit) is identical to the negative horizontal 
gradient of the moon. This rotational gradient is well within 
the dynamic range of the sensor and together with the horizontal 
gravity gradient of the moon will be part of the data set that 
will be used to produce not only the gravity field, but also the 
S/C orientation and motion (See Appendix D) . 
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(b) The gradients of the S/C attitude control system 
rotational rates will also be sensed. If momentum wheels are 
used and the rate errors are 2 microradians per second, the 
resulting error signal will be 0.004 Eotvos. 

This error is negligible for a sensor with a sensitivity of 
1 Eotvos. The uncertainties in the attitude control rate errors 
would have to be greater than 30 microradians per second to 
generate an error signal egual to the GEMS sensitivity. 

5.4 Attitude Control Accelerations 

The GEMS does not have first order sensitivity to S/C 
rotational or linear accelerations because of the common-mode 
rejection inherent in the design. 

5.5 Sensor Thermal Control 

5.5.1 Objective 

A thermal control study was conducted to determine the 
feasibility of maintaining RGG case temperature time gradient 
(rate of change) within +0.06 C/hr for a range of assumed S/C 
environmental temperatures, and temperature time gradients. 

5.5.2 Assumptions 

The assumed S/C environmental temperature and time gradient 
ranges for operation of the RGG within specification were 0 to 
55 C and 3 C/hr to 20 C/hr, respectively. In addition to the 
above assumptions, the orbit time was assumed to be 105 minutes 
and the worst case of 12 W internal RGG power dissipation was 

used. 


5.5.3 Preliminary Analysis 


A first look at the RGG temperature-control heater-power 
requirements indicated that the worst case boundary conditions 
ares (a) maximum ambient temperature range, (b) maximum 
environmental temperature rate of change, and (c) an orbit with a 
S/C sunlit-to-shadow ratio of one. 


The RGG nominal operating (case) temperature can be shifted 
periodically to maintain it at or slightly above the peak S/C 
ambient temperature. It is feasible to do this by program or 
command control. Thus, long term shifts in S/C ambient 
temperature as might be caused by shifts in the sun/shadow ratio 
or failure in S/C thermal control systems would not significantly 
increase RGG heater power requirements. 
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5.5.4 Thermal Design Approach 


The required thermal control of the RGG can be achieved by 
utilizing conventional thermal control techniques. Conduction 
isolation (conductivity less than 0.5) will be employed to 
minimize conduction effects of the S/C environment. High 
emittance thermal finish (emissivity greater than 0.85) is 
required to reject the 12 w internal dissioation of the unit 
while maintaining case temperatures below 61 C. Thess passive 
techniques will be utilized in conjunction with proportionally 
controlled heaters to null out case transient temperature 
variations. 

The combination of 12 watts internal dissipation and 
+.06 C/hour temperature stability requirement greatly influenced 
the thermal design of the unit. The present 0.18 square meter 
radiative surface area has to be approximately doubled to prevent 
the RGG bulk temperature from exceeding 61 C when exposed to the 
S/C maximum environmental temperature of 55 C. The required 
increase in radiative surface area may be accomplished by 
addition of fins to the ends of the RGG unit. 

For the final analysis, four maximum S/C environmental 
temperatures were selected: 55, 45, 30 and 15 C. The first 
represents the worst case maximum temperature at which the RGG 
must operate within specification. This is also the upper limit 
for the RGG, because of the temperature rist within the unit and 
the temperature limits of the Mil-Spec electronic components. 
The last is the lowest anticipated orbit-maximum S/C 
environmental operating temperature. The others are convenient 
intermediate values. Three values of temperature time gradient 
were used in the analysis: 3, 11, and 20 C/hour. Other factors 
used have been stated above. It should be emphasized that these 
include the worst case values for RGG internal power dissipation 
and sunlit-to-shadow ratio. 

5.5.5 Results 


The results of the thermal analysis are listed in Table 1. 
The first two columns are the selected parametric values for 
temperature rate of change and maximum environmental temperature. 
The third column contains the calculated values of minimum 
environmental temperature. The fourth column in Table 1 shows the 
calculated values of RGG bulk temperature, and the last two 
columns list the maximum and average heater power required for 
each combination of parametric conditions. 

The results are plotted in Figures 26 and 27. Figure 26 
shows the variation in S/C internal environmental temperature 
over the orbit period for each combination of parametric values. 
Figure 27 shows the variations in required average and peak 



1 

\ 



Table 1. RGG Thermal Analysis Results 
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ldT/A0| E 

(°C/hr) 

ENVIRONMENTAL 
TEMP, °C 

RGG BULK 
TEMP, °C 

HEATER POWER, W 

MAX 

MIN 

MAX 

AV6 

3 

15 

12.4 

22.3 

4.0 

2.0 


30 

27.4 

37.8 

5.0 

2.5 


45 

42.4 

51.0 

5.3 

2.7 


56 

52.4 

60.2 

5.6 

2.8 

11 

15 

5.4 

22.3 

13.8 

6.9 


30 

20.4 

37.8 

17.5 

8.8 


46 

35.4 

51.0 

18.8 

9.4 


66 

45.4 

60.2 

20.3 

10.2 

20 

16 

-2.5 

22.3 

23.9 

12.0 


30 

12.5 

37.8 

30.2 

15.1 


46 

27.6 

51.0 

32.6 

16.3 


66 

37.5 

60.2 

35.2 

17.6 


NOTES: 

1. lATMfllg- RATE OF CHANGE OF ENVIRONMENTAL TEMPERATURE 

2. ORBITAL PERIOD - 106 min 


3. RGO RADIATIVE SURFACE AREA - 0.36 m 2 

4. RGG INTERNAL POWER DISSIPATION - 12.0 W 
6. CONDUCTANCE FROM RGG TO S/C - 041 W/°C 


heater power over the range of maximum S/C internal 
environmental temperatures selected and for each selected value 
of temperature rate of change. Corresponding RGG bulk 
temperatures can also be read from the same curves using the 
upper abscissa scale. The peak power requirement occurs at the 
minimum S/C temperature and heater power decreases to zero at the 
maximum S/C temperature. The average power required is one half 
of the peak power. 
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Figure 26. Spacecraft Internal Environmental 
Temperature vs Orbit Time 
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Figure 27. RGG Bulk Temperature and Heater Power vs 
Maximum S/C Internal Environmental Temp 
and Temperature Change Rate 
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From the table and the curves it is seen that most 
significant factor dictating the heater power required to null 
out bulk case temperature variations is the S/C environment 
temperature change rate. Average (over one orbit) heater power 
consumption ranges from 2 to 2.8 W for a 3 C/hr change rate, 6.9 
to 10.2 W for an 11 C/hr rate and 12.0 to 17.6 W for a 20 C/hour 
rate. 


5.5.6 Conclusions 


Thermal control of the RGG to maintain case temperature 
within +0.06 C is feasible with a radiation dominated design, 
utilizing proportional controlled heaters to null out S/C 
internal envir*. .oriental temperature changes. Location of the RGG 
in a region not likely to have high environmental temperature 
change rates is recommended to minimize heater power consumption. 

5.6 Sensor Dynamic Characteristics 

The RGG instrument has a significant amount of angular 
momentum (30 million gm-cm squared/second), but due to the speed 
control of one part in 10 million, this should not disturb the 
S/C attitude. A brief consideration of the momentum effect on 
attitude change control indicates that there will be no increase 
in thrust requirement, if the RGG momentum is taken into account 
in the attitude control equations and mechanization. Further 
study is required to fully evaluate momentum effects. 
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6.0 INTERFACE DEFINITIONS 

6.1 M^chanical/Weiqht/Volume 

(a) Sensor Mechanical : The mechanical configuration of 

the RGG sensor is shown in Figure 2. While the exact interface 
will be a function of the S/C, the unit can be mounted using the 
existing center flange. The addition of radiation fins as 

discussed in Section 5.5 will increase the length slightly. 
Nominal S/C alignment accuracy is adequate (see Section 5.2). 

(b) Electronics Mechanical : The configuration of the 

external electronics package for one RGG sensor is shown in the 
layout and mounting drawing of Figure 25. The case, except for 
the cover, and the structure for supporting the electronic 
subassemblies are machined from a single aluminum casting. This 
produces a very rugged unit and also serves as a good heat 
conductor between the electronics and the S/C structure. 

(c) Sensor Weight : The RGG sensor including its 

internal electronics weighs approximately 16 Kg. 

(d) Electronics Weight : The external electronics for 

one RGG sensor will weigh about 3 Kg. 

(e) Sensor Volume : The sensor is approximately 30 cm 

long by 23 cm in diameter. 

(f) Electronics Volume : The electronics unit is 

37x15x8 cm. 


6.2 Thermal 

(a) Sensor : The RGG will have its own thermal system 
to control temperature to within +.06 C .of the nominal operating 
temperature set by command signal. This control will be 
effective over a spacecraft environmental temperature range from 
below 0 C to 55 C, and for temperature rates of change up to 
20 C/hr. The cost of wide and rapid temperature changes, 
however, is the requirement for increased heater power to 
maintain the RGG case temperature within the +0.06 C tolerance. 
This power requirement will be kept within bounds by occasionally 
changing the nominal RGG temperature operating point 
automatically or by command as the S/C interior environmental 
temperature changes due to changes in the sun/shade ratio or 
other factors. Spin speed must be changed when the nominal 
temperature operating point is changed. The heater power 
requirement can be reduced to a very nominal value of 2 to 3 
watts average by locating the RGG in an area of the S/C which has 
a low temperature rate of change. The sensor-to-spacecraf t 
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thermal conductivity will be approximately 0.5 W/C. The thermal 
study results are discussed in more detail in Section 5.5. 

(b) Electronics : The Hughes standard space-qualified 

electronics design and fabrication techniques allow operation 
within specification over a range of -7 C to 55 C and survival 
without damage over a range of -40 C to 71 C. Thermal 
conductivity between the external electronics unit and the S/C is 
made as low as practicable. 

(c) Thermal Power : 

RGG: 6 to 12 W (Probably near 6 W) 

External Electronics: 6 W 

Heater: 2 to 18 w avg (See (a) above and 

Section 5.5) 


6.3 Power 


The GEMS for space applications operates from raw 28vdc 
power with power switching provided by the S/C. Total 
consumption is between 14 and 36 W average over one orbit, or 
between 16 and 54 W peak at the coldest period in the orbit. 
Both ranges are calculated using worst case conditions. The 
actual power in each case will depend on changes expected to be 
made in the spin motors (-6 W) and primarily on location of the 
RGG in the S/C, effecting temperature rate of change and thus 
heater power (See Section 6.2 (a) above and Section 5.5). No 
pulsed power modes are used. 

6.4 Commund 


Command Requirements for GEMS are listed in Table 2. 

6.5 Data 


Following are data register and data rate requirements for 

GEMS. 


(a) Science Data : Table 3. The output science data 
should be identified with the lunar surface position over which 
it was gathered to an accuracy equal to 0.1 to 0.3 times the S/C 
nominal altitude above the lunar surface. 

(b) RGG Engineering Data : Table 4. 


(c) External Electronics Engineering Data i Table 5. 
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Table 2. Command Words 
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ITEM 

NO. 

COMMAND DATA ITEM 

WORDS 

BITS 

1 

MOTOR SPEED (PRESET OCTANT COUNT) 

1 

17 

2 

MOTOR SPEED VARIATION LIMIT 

1 

8 

3 

MOTOR STOP 

1 

1 

4 

MOTOR DRIVE SERVO VELOCITY GAIN 

1 

10 

5 

MOTOR DRIVE SERVO PHASE GAIN 

1 

10 

6 

NOMINAL RGG OPERATING TEMPERATURE 

1 

7 

7 

RGG TEMPERATURE TOLERANCE LIMIT 

1 

4 

8 

RGG DATA READ 

1 

2 

9 

CAPACITOR DRIVER 

1 

1 

10 

RGG MODE CONTROL 

1 

4 

11 

ENGINEERING DATA RATE 

1 

1 


TOTAL 

11 

65 


Table 3. Science Data 
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ITEM 

NO. 

SCIENCE DATA ITEM 

WORDS 

BITS 

PER WORD 

DATA RATE j 

WORDS/MI N 

BITS/MIN 

1 

GRAVITY GRADIENT 
MEASUREMENTS: 

EIGHT “OCTANT COUNTS" OF 
10 MHz REFERENCE CLOCK 
(MAY INCLUDE TEST 
REFERENCE TORQUE) 

8 

8x20^ 

8 x 16 

160 x 16 


TOTAL 

S 

1S2 

(bit!) 

128 

2660 


m 


TWENTY BITS INCLUDES 19 BITS FOR GRAVITY GRADIENT DATA AND 1 BIT FOR 
WIDE SLOT MARKER (PHASE REFERENCE). 




Table 4. 


RGG Engineering Data 
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ITEM 

NO, 

Dee cmnucCDUUn 


BITS 

PER WORD 

DATA RATE j 

Kuu tNulWCtnrniu 

DATA ITEM 

WORDS 

WORDS/MIN 

BITS/MIN 

1 

UNREGULATED VOLTAGE 

1 

8 

1/12 121 

8/96 

2 

REGULATED VOLTAGE 

1 

8 

1/12 

8/96 

3 

RF VOLTAGE 

1 

8 

1/12 

8/96 

4 

TEMPERATURE SENSOR NO. 1 

1 

10 

1/12 

10/120 

S 

TEMPERATURE SENSOR NO.2 

1 

10 

1/12 

10/120 

b 

REFERENCE TWO-OMEGA 
SQUARE WAVE VOLTAGE 

1 

8 

1/12 

8/96 


TOTAL 

6 

52 (bits) 

6/72 

52/624 | 


l2] THE NORMAL RATE OF 1 WORD PER MIN IS MORE THAN ADEQUATE, EXCEPT WHEN 
MONITORING CHANGES IN OPERATING PARAMETERS DURING AND IMMEDIATELY 
AFTER OPERATIONAL COMMAND INPUTS. THEN A RATE OF AT LEAST 12 WORDS 
PER MIN IS DESIRABLE. 



Table 5. External Electronics Engineering Data 
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ITEM 

EXTERNAL ELECTRONICS 

WORDS 

BITS 

DATA RATE | 

NO. 

ENGINEERING DATA ITEM 

PER WORD 

WORDS/MIN 

BITS/MIN 

1 

REFERENCE 2 w 
SQUARE WAVE VOLTAGE 

1 

8 

1/12 131 

8/96 

2 

FOUR POWER SUPPLY VOLTAGES 

4 

8 

4x (1/12) 

4 x(8/96) 

3 

MOTOR CONTROL FIELD 

1 

8 

1/12 

8/96 

4 

MOTOR REFERENCE FIELD 

1 

8 

1/12 

8/96 

5 

MOTOR SPEED ERROR 

1 

1 

1/12 

1/12 

6 

MOTOR STOP 

1 

1 

1/12 

1/12 

7 

WIDE SLOT MARKER FAIL 

1 

1 

1/12 

1/12 

8 

VELOCITY ERROR 
(8 QUADRANTS) 

8 

5 

8 x (1/12) 

8 x (5/60) 

9 

FM SIGNAL LEVEL 

1 

8 

1/12 

8/96 

10 

Ki DATA 

8 

12 

8 x (1/12) 

8x02/144) 

11 

STATOR TEMPERATURE 

1 

10 

1/12 

10/120 

12 

AMBIENT TEMPERATURE 

1 

10 

1/12 

10/120 

13 

STATOR TEMPERATURE ERROR 

1 

1 

1/12 

1/12 

14 

HEATER CURRENT 

1 

8 

1/12 

8/96 

15 

RGG LOGIC CODE 

1 

4 

1/12 

4/48 

16 

PHOTO PICKOFF SIGNAL 

1 

8 

1/12 

8/96 


TOTAL 

33 

244 (bits) 

33/396 

244/2928 


131 (SEE FOOTNOTE OF TABLE 3) 
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7.0 PERFORMANCE PARAMETERS 


7.1 Sensitivity 


The GEMS has a design sensitivity of 1 Eotvos for a 10 
second integration period. 

7.2 Dynamic Range 

The GEMS is capable of operating over a gravity gradient 
range of +50,000 Eotvos to -50,000 Eotvos. 

7.3 Quantizatio n 

By using the highly accurate digital speed control servo in 
each RGG to demodulate and digitize the analog sensor output, the 
GEMS converts the gravity gradient signal to a digital word with 
a precision of 0.2 Eotvos (19 bits). 

7.4 Alignment 

The RGG does not need to be accurately aligned, because 
ground data processing of GEMS data results in a knowledge of 
alignment with respect to local vertical of + 0.001 radian. 

7.5 Spin Speed 

The nominal spin speed of the RGG is 1050 rpm (17.5 rps) 
with the capability to control the speed at a commanded value 
from 500 rpm (0 rpm with design modification) to 1500 rpm. 

7.6 Spin Speed Accuracy 

The RGG spin speed is controlled to an accuracy of one part 
in 10 million. The short term noise is less than 20 microradians 
and the long term phase stability is better than one part in 1 
billion. See Section 3.3.2 for details on the RGG spin speed 
control. 


7,7 Temperature 

RGG bulk (case) temperature is maintained within +0.06 C 
about a commanded case operating temperature which is set ”at or 
slightly above the sensed S/C maximum internal environmental 
temperature at the location of the RGG. See Section 5.5 for 
details on the RGG thermal control. 
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7.8 Sensitivity to Spurious External Fields 


The RGG is fabricated from non-magnetic materials to avoid 
magnetic field error sensitivity. In the original design, an 
inadvertent closed loop conductive path through each set of arm 
plates and end masses has resulted in a higher magnetic field 
error sensitivity than was expected. RGG prototype No. 1 has a 
magnetic field error sensitivity of approximately 20,000 Eotvos 
per Gauss squared. A design modification has already been 
incorporated in RGG No. 2 to break the closed conductive path. 
It is anticipated that this will reduce the magnetic field error 
sensitivity by the required amount. There will still remain an 
eddy current effect which j difficult to calculate. If the 
sensitivity is still above the level of 1000 Eotvos per Gauss 
squared, and if the sensor could be subject to magnetic field 
variations greater than 0.01 Gauss from other equipment in the 
spacecraft, magnetic shielding will be used. This has been 
proven in the laboratory to be completely successful in reducing 
the error to acceptable limits, even with RGG No. l's high 
magnetic field error sensitivity and operating in the earth's 
magnetic field. Thus, magnetic field error sensitivity will not 
be a problem in the LPO application. 

The RGG is not sensitive to electrostatic fields nor to 
nuclear radiation. Sensitivity to mass shifts in the S/C can be 
minimized by locating the RGG away from large moving masses such 
as platforms, neglecting the gravity gradient data during periods 
of large mass shifts, if infrequent or short in duration, or by 
compensating for the shifts in the data if needed. 

The external electronics has the same degree of sensitivity 
to nuclear radiation as conventional space-qualified analog and 
digital circuitry. Hardened devices and design techniques can be 
used if required. 

7.9 Generation of Interference Fields 


The RGG does not utilize nor produce nuclear radiation. It 
produces only the normal electromagnetic fields produced by 
current flow in cables external to the RGG housing and 
electronics unit enclosure. These can be minimized by proper 
design and cable shielding. 

The radio frequency (RP) link is between the RGG rotor and 
stator and thus is totally inclosed within the RGG metal housing 
(stator) which does not leak RF. Therefore, the RF link does not 
produce radio frequency interference (RFI) external to the RGG. 
Even within the RGG the RF link produces no significant effect on 
the internal electronics circuitry including the low level 
preamplifier for the sensor signal from the piezoelectric 
transducer. 
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8.0 


SAFETY 


There are no high voltages, radioactive materials or other 
hazardous items in GEMS. 

The RGG is designed to withstand up to 50 g shock loads in 
the radial direction (perpendicular to the spin axis) and up to 
100 g in the longitudinal direction (parallel with the spin axis) 
under non-operating conditions (for ground handling and launch). 
These factors can be increased, if necessary, by making a minor 
design change. 

No problems are anticipated in operating the GEMS in the 
very benign space environment as it was designed for operation in 
aircraft and other vehicles with more severe environments. 

9.0 MISSION 


(TBD) 



APPENDIX C 


LPOSIM FORTRAN PROGRAM 


SIMULATES A LUNAR POLAR ORBITER GRAVITY GRADIENT EXPERIMENT 




c 

C LPOSIM. FOR 
C 

C SIMULATES LUNAR POLAR ORBITER GRAVITY GRADIENT EXPERIMENT 
C 

C 30 JUNE 1976 VERSION 
C 

C DR. ROBERT L. FORWARD, SENIOR SCIENTIST 

C HUGHES RESEARCH LABS, 3011 MALIBU CANYON RD. , MALIBU, CA 90265 
C (213)456-6411 
C 

C CALCULATES THE VARIOUS GRAVITY GRADIENTS AND COMBINATIONS OF 
C GRAVITY GRADIENTS TO BE EXPECTED ON THE LPO. THE LPO 13 IN A 
C SIMPLE ELLIPTICAL ORBIT ABOUT A SPHERICAL, INERT I ALLY ROTATING, 
C BUT NON-TRANSLATING MOON. 

C THE MOON HAS VARIOUS MASCONS ON IT THAT ARE MOVED UNDER THE LPO 
C ORBIT AS THE MOON ROTATES. 

C 

C DIMENSION THE MASS POINT SET AND OPEN DATA STORAGE FILES. 

C 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DIMENSION AM (20) ,XM(20) ,YM(20) , ZM ( 20 ) 

C 

OPEN (UNIT=20,FILE='GVV. DAT') 

OPEN (UNIT=21,FILE='GAA. DAT') 

OPEN (UNIT=22,FILE='GCC. DAT') 

OPEN (UNIT=23 , FILE= 'GVA.DAT ' ) 

OPEN (UNIT=24,FILE='GAC. DAT') 

OPEN (UNIT*25 , DEVICE* ' DSK ' , FILE* ' GCV . DAT ' ) 

OPEN (UNIT* 2 6 , DEVICE* 'DSK ' , FILE= ' GVMA . DAT ' ) 

OPEN (UNIT=27 , DEVICE* ’ DSK ' , FILE* ' GAMC . DAT ' ) 

OPEN (UNIT*28, DEVICE*’ DSK ', FILE* 'GCMV.DAT') 

OPEN (UNIT*29,DEVICE*'DSK' , FILE* 'GTRA.DAT' ) 

OPEN (UNIT*30, DEVICE*' DSK ’, FILE* 'ALT. DAT') 

OPEN (UN IT* 3 1 , DEVICE* ' DSK ’ , FI LE* ' TI ME . DAT ' ) 

C 

C SET UP CONSTANTS. 

C G * GRAVITATIONAL CONSTANT IN KM*3/KG*SEC*2 

C A * LUNAR RADIUS IN KM 

C AM , X , Y , Z * MASS AND POSITION OF MOON 

C WL* LUNAR SIDEREAL ROTATION IN RAD/SEC 

C 

PI*3. 1415926 

G-6.67E-20 

A-1737.9 

AM(1) *7 . 349E+22 

XM(1) *0 

YM(1) *0 

ZM (1) *0 

WL-2.6618E-6 
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SET UP ORBITAL PARAMETERS FOR SPACECRAFT 
E * ORBITAL ECCENTRICITY 

H = MINIMUM ORBITAL ALTITUDE IN KM (PERILUNE) 

P = SPACECRAFT ANGULAR MOMENTUM IN KG*KM* 2/SEC 
RP = RADIUS CONSTANT IN ORBITAL EQUATION IN KM 
TH0= INITIAL ANGLE FOR ORBIT IN RAD 

E=0. 0000 
TH0=-0.25*PI 
H=60 . 

RP= (A+H) * (1+E) 

P=SQRT(G*AM(1) *RP) 

START SPACECRAFT ORBITING MOON 

TH = SPACECRAFT ANGLE ALONG ORBIT IN RAD 
DT = TIME INCREMENT IN SEC 

RS = SPACECRAFT RADIUS FROM LUNAR CENTER IN KM 
ALT® SPACECRAFT ALTITUDE IN KM 

WS = SPACECRAFT INERTIAL ANGULAR RATE IN RmD/SEC 
XS , YS , ZS ■ SPACECRAFT POSITION IN KM 

TH=TH0 

DT=10 

DO 10 T=-24000, 24000, DT 
RS=RP/(l+E*COS(TH) ) 

ALT-RS-A 

WS»P/(RS*RS) 

XS«RS*COS(TH) 

YS»0 

ZS-RS*SIN(TH) 

GENERATE MASS POINT SET 

AM (I) * MASS OF Ith MASS POINT 

XM(I) ,YM(I) ,ZM(I)« POSITION OF MASS POINT IN INERTIAL FRAME 

SPACE THEM 2 DEG « 0.035 RAD APART 
CD-COS (0.035) 

SD*SIN (0.035) 

AND MOVE THEM AROUND WITH THE LUNAR ROTATION. 

CWLT-COS (WL*T) 

SWLT«SIN(WL*T) 
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MASCON DIRECTLY ALONG X AXIS 
AM (2) =5 . OE+16 
XM (2) «A*CWLT 
YM(2) =A*SWLT 
ZM (2) *0 

C MASCON ON EQUATOR-WEST OF #2 
AM ( 3) =0 

XM ( 3) =A*CD*CWLT+A*SD*SWLT 
YM ( 3) =-A*SD*CWLT+A*CD*SWLT 
ZM ( 3) *0 

C MASCON BELOW #2 
AM (4 ) =5 . GE+16 
XM ( 4 ) =A*CD*CWLT 
YM(4)=A*CD*SWLT 
ZM ( 4 ) =-A*SD 

C MASCON ALONG -X AXIS (FARSIDE) 
AM(5)=5. OE+16 
XM ( 5) =-A*CWLT 
YM ( 5) =-A*SWLT 
ZM ( 5) =0 

C MASCON ON EQUATOR-WEST OF #4 
AM (6) =5 . OE+16 

XM ( 6) =-A*CD*CWLT-A*SD*SWLT 
YM(6)=-A*CD*SWLT+A*SD*CWLT 
ZM (6 ) =0 

C MASCON BELOW #4 
AM ( 7 ) =0 

XM(7)=-A*CD*CWLT 
YM ( 7 ) »-A*CD*SWLT 
ZM (7) =A*SD 

C MASCON LUNAR NORTH POLE 
AM(8)=5. OE+16 
XM(8) *0 
YM (8) =0 
ZM (8) =A 

C MASCON SOUTH OF #8 TOWARD #2 
AM (9) *0 

XM(9)»A*SD*CWLT 

YM(9)»A*SD*SWLT 

ZM(9)«A*CD 
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c 

C CALCULATE TOTAL GRAVITY GRADIENT OF ENTIRE MASS POINT SET 
C IN INERTIAL REFERENCE FRAME. 

C 

GXX=0 

GYY=0 

GZZ-0 

GXY=0 

GYZ-0 

GZX=0 

DO 20 1-1,9 
GM-G*AM(I) 

X=XM ( I ) -XS 
Y-YM ( I ) -YS 
Z=ZM ( I ) -ZS 
R-SQRT (X*X+Y*Y+Z*Z) 

GRR-GM/ ( R*R*R) 

GXX-GXX+GRR* (-1+3*X*X/(R*R) ) 

GYY=GYY+GRR*(-1+3*Y*Y/(R*R) ) 

GZZ-GZZ+GRR* (2-3* (X*X+Y*Y) / (R*R) ) 

GXY*GXY+GRR*3*X*Y/ (R*R) 

GYZ=GYZ-GRR*3*Y*Z/(R*R) 

GZX=GZX-GRR*3*Z*X/(R*R) 

20 CONTINUE 
C 

C ROTATE TOTAL GRADIENT CALCULATED IN INERTIAL REFERENCE FRAME TO 
C A LUNAR VERTICAL REFERENCE FRAME, V * VERTICAL, A * ALONG TRACK 
C AND C = CROSS TRACK. 

C 

CWCT-COS (TH) 

SWCT-SIN (TH) 

GVV«GXX*CWCT*CWCT“2*GZX*SWCT*CWCT+GZZ*SWCT*3WCT 
G AA-GX X * S WCT * SWCT+ 2 *G Z X * SWCT*CWCT+G Z Z * CWCT * CWCT 
GCC-GYY 

GVA- (GXX-GZZ) *SWCT*CWCT+GZX* (CWCT*CWCT-SWCT*SWCT) 

GAC— GXY*SWCT-GYZ*CWCT 
GCV— GX Y *CWCT+G Y Z * SWCT 
C 

C NOW PUT IN EFFECTIVE GRADIENTS OF THE SPACECRAFT ANGULAR RATES 
C 

WV-0 

WA-0 

WC-WS 

GW«GW+WA*WA+WC*WC 
GAA«GAA+WV*WV+WC*WC 
GCC-GCC+WV*WV+WA *WA 
GVA«GVA+WV*WA 
GAC-GAC+WA*WC 
GCV-GCV+WC * WV 
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C 

C PUT IN ANGULAR PERTURBATIONS ABOUT THE LUNAR VERTICAL 
C COORDINATES (USE SMALL ANGLES ONLY SINCE THE ROTATION 
C MATRICES USED ARF, NOT INDEPENDENT) 

C 

C PUT IN SMALL PERTURBATION ABOUT VERTICAL AXIS 
C VP = PERTURBATION ANGLE IN RADIANS (SMALL) 

C 

VP=0 

CV^COS (VP) 

SV=SIN (VP) 

PVV=GVV 

PAA=GAA*CV*CV+2*GAC*SV*CV+GCC*SV*SV 
PCC-GCC*CV*CV-2*GAC*SV*CV L GAA*SV*SV 
PVA=GVA*CV+GCV*SV 

PAC=GAC* (CV*CV-SV*SV) + (GCC-GAA) *SV*SV 
PCV=GCV*CV-GVA*SV 

PUT IN SMALL PERTURBATION ABOUT ALONG TRACK AXIS 
AP=PERTURBATION ANGLE IN RADIANS (SMALL) 


AP=0 

CA=COS (AP) 

SA=SIN (AP) 

GVV=PVV*CA*CA+2*PCV*SA*CA+PCC*SA*SA 

GAA=PAA 

GCC-PCC*CA*CA-2*PCV*SA*CA+PVV*SA*SA 

GVA=PVA*CA+PAC*SA 

GAC=PAC*CA-PVA*SA 

GCV=PCV* (CA*CA-SA*SA) + (PCC-PVV) *SA*CA 

PUT IN SMALL PERTURBATION ABOUT CROSS TRACK AXIS 
CP=PERTURBATION ANGLE IN RADIANS (SMALL) 


CP=0 

CC=COS(CP) 

SC=SIN (CP) 

PVV=GVV*CC*CC+2*GVA*SC*CC+GAA*SC*SC 

PAA=GAA*CC*CC-2*GVA*SC*CC+GVV*SC*SC 

PCC=GCC 

PVA=GVA* (CC*CC-SC*SC) + (GAA-GVV) *SC*CC 

PAC=GAC*CC-GCV*SC 

PCV=GCV*CC+GAC*SC 
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c 

C CONVERT GRADIENTS TO EOTVOS UNITS. 

C 

GVV=PVV*1 . oe+9 
GAA=PAA*1 . OE + 9 
GCC=PCC*1 .OE+9 
GVA=PVA*1 .OE+9 
GAC=PAC*1 .OE+9 
GCV=PCV*1 .OE+9 
C 

C CALCULATE IN-PHASE OUTPUTS OF ROTATING GRAVITY G RADIOMETERS 
C AND TRACE OF GRADIENT TENSOR. 

C (TRACE SHOULD EQUAL TWICE THE SQUARE OF THE ANGULAR RATES.) 
C 

GVMA= (GVV-GAA) /2 . 0 
GAMC= (GAA-GCC ) /2 . 0 
GCMV= (GCC-GVV) /2 . 0 
GTRA=GVV+GAA+GCC 
C 

C WRITE RESULTS IN FILES 
C 

WRITE (20,1) GW 
WRITE (21,1) GAA 
WRITE (22,1) GCC 
WRITE (23,1) GVA 
WRITE (24,1) GAC 
WRITE (25, 1)GCV 
WRITE (26 ,1)GVMA 
WRITE (27, 1)GAMC 
WRITE ( 28, 1)GCMV 
WRITE (29, 1)GTRA 
WRITE (30,1) ALT 
WRITE ( 3 1 , 1 ) T 
1 FORMAT ( F10 . 2) 

C 

C CONTINUE MOTION ALONG ORBIT 
C 

TH=TH+WS*DT 
10 CONTINUE 
C 

C STOP WHEN FINISHED 
C 

STOP 

END 
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LUNAR POLAR ORBITER GRAVITY GRADIENT EXPERIMENT 
SIMULATION PLOTS 


To examine the behavior o£ a gravity gradient measurement 
system on board an LPO spacecraft, we used the simulation program 
in Appendix C. 

In our first run with the program, we assumed a circular 
60 Km orbit for the spacecraft, and did not introduce any pertur- 
bations in the spacecraft attitude. We inserted three mascons. 
One mascon was on the near side, one at the north pole, and one 
on the far side. The plots from this simulation are shown in the 
first two figures. The six outputs of the sensors (see Fig- 
ure {1} and {2}) were flat, except where they went across the 
mascons. in Figure {1} we see that all the "principal gradient 
difference" outputs of the RGG sensor o have significant biases. 
The bias level on Gvma or (Gvv-Gaa)/2 is 1266 Eotvos, on Game 
or (Gaa-Gcc)/2 is 422 Eotvos, and on Gcmv or (Gcc-Gvv)/2 is 
-1688 Eotvos. The three cross gradient outputs (See Figure {2}) 
were zero except when they crossed a mascon. 











From these sensor outputs we can calculate the three princi- 
pal gradient tensor components (See Figure {3}) and the trace of 
the gradient tensor. 

For the zero eccentricity, zero attitude error case, the 
principal vertical gradient Gw bias is 2532 Eotvos, the princi- 
pal cross track gradient Gcc bias is -844 Eotvos and the princi- 
pal along track gradient Gaa bias is zero. This last occurs 
because the along track gravity gradient is exactly canceled by 
the rotational rate gradient of the spacecraft rotating once per 
orbit to maintain a local vertical orientation. 

Because the spacecraft is vertically stabilized, it is 
rotating with respect to inertial space. Since the trace of the 
gradient tensor Gtra is independent of the lunar or mascon 
gravity fields and is only equal to the square of the angular 
rate of the spacecraft, it was constant at 1688 Eotvos. 





We then introduced an eccentricity of e®0.005 into the 
orbit. This caused the spacecraft altitude to vary from 60 Km on 
the near side to 78 Km on the far side. With this altitude vari- 
ation we get variations in the measured gradient of tens of Eot- 
vos because of the background gravity gradient of the moon. But 
as can be seen in Figure {4} through Figure {6}, these variations 
occur very slowly compared to the more rapid changes in the grad- 
ient as we cross over a mascon and their strict once per orbit 
periodicity should allow us to remove them from the sensor data. 

Note, however, that not all the gradiometer outputs are 
affected. As we can see in Figure {5}, the sensor outputs that 
give a measure of the cross gradients do not have any shift in 
their bias level due to the eccentricity of the orbit since their 
initial bias was zero. Thus, even for an elliptical orbit, we 
can use this data just as it comes from the sensor without any 
preprocessing and integrate it to get the gravity. 







TIME, mc 


Figure (6) - Derived Principal Gradients 
60 Km Pcrilune (e * 0.005) 








If we have three sensors and can obtain the trace of the 
gradient tensor , we find that the trace has the same variation as 
the altitude (See Figure {7}), but it has no gravity information 
in it since the trace of the gravity part of the gradient tensor 
is zero. The variation seen in the trace of the gradient tensor 
is just the angular rate variation due to the ellipticity of the 
orbit and is a direct measure of the ellipticity and altitude. 
The amplitude of the variations in the trace are about 4 Eotvos 
for a 1 Km change in altitude. Since we can curve fit over many 
minutes (say 1000 sec) , then the accuracy of the curve fit would 
be 0.1 Eotvos and we see that we can use the output of the trace 
of the gravity gradient measurement system to estimate the space- 
craft altitude variations to better than 25 m. 


} 
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Figure {7} - Spacecraft Altitude and Gradient Tensor Trace 
60 Km Perilune (e * 0.005) 
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It is instructive to examine the very low frequency harmonic 
portions of the many gradiometer outputs, since they have signi- 
ficantly different responses to the altitude variation and angu- 
lar rate variations that are caused by the ellipticity of the 
orbit. 

In a simple elliptical two body orbit, the radius R of the 
orbit varies as 


1/R = (1 + e cos (a)) /A 

where A is the average orbital radius 

e is the eccentricity of the orbit 
a is the angle along the orbit 

The square of the orbital angular rate Wcc of the spacecraft 
for this simple orbit varies as: 

Wcc = Wc x Wc = Grr (1 + e cos (a)) 

3 

where Grr = GM/R is the gravity gradient amplitude. 

Wc * orbital angular rate 

Since the gravity gradient amplitude Grr varies as the 
inverse cube of R , if the eccentricity of the orbit is small, 
this can be approximated by: 

Grr = Ggg (1 + 3e cos (a)) 


3 

where Ggg * GM/A is the average gravity gradient amplitude 

In the same manner, the square of the orbital angular rate 
has the variation: 


Wcc * Ggg (1 + 4e cos (a)) 

Because the LPO spacecraft is stabilized to the Lunar verti- 
cal by a control system sensing the Lunar horizon, the spacecraft 
is rotating with respect to the inertial space at an instantane- 
ous rate that is identical to the instantaneous orbital angular 
rate. The gradiometer instruments will measure the gradients of 
both the Lunar gravity field (which will vary with the change in 
orbital radius) , and the gradient of the angular rate (which var- 
ies in a different manner with the motion along the orbit) . 
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After a bit of calculation it can be shown that the low fre- 
quency components of the six gradiometer outputs have the varia- 
tion : 


Gvma = 1.5 Grr = Ggg 
Game = 0.5 Wcc = Ggg 
Gcmv =-1.5 Grr -0.5 Wcc 
Gva = Gac 


(1.5 + 4 . 5e cos (a) ) 

(0.5 + 2 . Oe cos (a) ) 

= Ggg (-2 -6.5e cos(a)) 
= Gcv = 0 


Notice that the Gvma gradiometer output does not contain 
any of the spaceerc c t angular rate in the cross track direction 
(the orbital axis direction) , thus its output is pure gravita- 
tional data. The low frequency variations seen in this output 
are then directly related to the orbital radius variation through 
the inverse cube law of the gravity gradient. 


From the gradiometer outputs we can derive other outputs 
such as the principal components of the gradient tensor: 


Gw = 

2Gr r 

+ Wcc = Ggg 

(3 

+10e 

cos (a) ) 

Gaa = 

-Grr 

+ Wcc = Ggg 

(0 

+ e 

cos (a) ) 

Gcc = 

-Grr 

= Ggg 

(1 

+ 3e 

cos (a) ) 


and the trace of the gradient tensor 

Gtra = 2 Wcc = Ggg (2 + 8e cos (a)) 

If we have a three axis gradiometer system, then the gradi- 
ent trace will be a very important output, for it will contain no 
gravity gradient signals at any harmonic order, so any variations 
seen in the trace - fast or slow - are due to angular rate varia- 
tions. Their presence in the trace will mean that similar varia- 
tions will be found in the other outputs of the system and that 
these output variations are due to angular rate variations and 
are not due to gravity anomalies. 
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We then repeated the simulation for a circular orbit at 
60 Km with the spacecraft attitude off about the along track 
direction (roll) by 1 degree (17 milliradians) . We now again get 
essentially flat traces except for the mascons (See Figures {8} 
through {10}). The RGG principal gradient difference outputs 
(See Figure {8}) are almost the same as they were for perfect 
attitude (the cosine of one degree is very close to one) and the 
cross gradient outputs (See Figure {9}) have zero bias except for 
Gcv , which now has a constant bias of -58 Eotvos due to the 
roll of one degree. Thus, by attributing this bias to a space- 
craft attitude error, we could use the cross gradient component 
bias to estimate the spacecraft attitude error to much better 
than a milliradian. 

Notice that in the plots of the cross gradient terms Gac 
and Gcv , the north pole mascon shows up as a slight signal 
whereas for perfect roll orientation (See Gac and Gcv in Fig- 
ure {2}), those cross gradient components remain zero during the 
polar mascon crossing. These errors could be eliminated in post 
data processing using the Gcv bias knowledge, but even if they 
are not, the effect is only to misplace the mascon by 600 meters 
for 60 Km altitude - a shift that is far less than the data reso- 
lution. 













I 
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Since it seemed from the plots and our previous experience 
that a great deal of information could be obtained from a single 
RGG sensor oriented along the orbit axis so that it measured the 
gradients in the orbital plane, we decided to generate some simu- 
lations using only the output of one RGG sensor. One of the 
questions was whether a single sensor oriented to sense the grad- 
ients along the orbital track would provide adequate resolution 
for mascons separated in the cross track direction. 

To look at this question we rearranged the orbital equations 
to give us traces of the gravity gradient output as the orbit 
passed over the front side of the moon in 20 successive orbits. 
On the front side we placed three mascons. One directly at the 
nearside point, one north of it (along the orbital track) and one 
to the west (in the cross track direction). We first spaced the 
mascons at 60 Km separation, the same distance as the altitude. 
The plots were then presented in a perspective view. 



We first needed something to compare the single sensor out- 
put with. In Figure 111} is plotted the vertical gradient of the 
vertical gravity Gw . Our previous simulations had shown that 
this gradient component was of high resolution and gave a good 
representation of the mass distribution. Examination of the plot 
shows that this component of the gravity gradient tensor can sep- 
arate the three mascons, both along track and cross track even 
when the spacecraft altitude is equal to the mascon spacing. 



Figure {11} - Vertical Gradient of Vertical Gravity 
Three Mascor with 60 Km Spacing 
60 Km Ciro iar Orbit 
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Now the vertical gradient of the vertical gravity Gw is 
not one of the normal outputs of a single RGG system, but has to 
be derived from the outputs of three gravity gradiometers. The 
nearest equivalent sensor output of a single RGG system is the 
sine output of the RGG which is the difference of the vertical 
gradient and the along track gradient (Gvv-Gaa)/2. This sensor 
output is plotted in Figure {12}. We see that it is a slightly 
distorted version of Gw . It has equivalent resolution in the 
cross track direction, and slightly better resolution in the 
along track direction. 
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Figure {12} - RGG Principal Gradient Difference Output 
Three Mascons at 60 Km Spacing 
60 Km Circular Orbit 
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Figure {14} is the integrated cross gradient output of a 
single RGG sensor. There is a nice rounded peak, but it is obvi- 
ous that the integrated cross gradient does not have the resolu- 
tion to see the three mascons when they are separated by a dis- 
tance equal to the spacecraft altitude. However, this integrated 
output gives an undistorted gravity contour map. Thus for opti- 
mum data processing, this undistorted but lower resolution inte- 
grated output should be used in conjunction with the higher reso- 
lution (but slightly distorted) principal gradient difference 
output to provide a high resolution, low distortion data set for 
determining underground mass distributions. 
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Figure {14} - Integrated RGG Cross Gradient Output 
Three Mascons at 60 Km Spacing 
60 Km Circular Orbit 
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We next wanted to study the effects of an elliptical orbit 
on the data resolution. We used the same three mascon set as in 
the previous plots, but now we spaced them 120 Km (4 degrees) 
apart. We chose an orbital eccentricity of e=G.005 (which gives 
an orbital altitude variation from 60 to 78 Km) and we placed the 
perilune of the ellipse right over the three mascon set. In Fig- 
ure {15} is shown the best that could be obtained from a three 
sensor system, the vertical gradient of the vertical gravity. 
Notice, that although there is a significant curvature to the 
data, the mascon peaks are still resolved. 



Figure {15} - Vertical Gradient of Vertical Gravity 
Three Mascons at 120 Km Spacing 
60 Km Perilune of e»0.005 Elliptical Orbit 
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A plot of the principal gradient difference output Gvma of 
a single RGG for the same orbital conditions is shown in Fig- 
ure {16}. It too has good sensitivity and resolution for the 
elliptical orbit case. 




Figure {16} - RGG Principal Gradient Difference Output 
Three Nascons at 120 Km Spacing 
60 Km Perilune of e-0.005 Elliptical Orbit 
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The other output of a single RGG system is the cross gradi- 
ent output, and as before, we can integrate it to obtain the 
vertical gravity. Although there might be some concern about the 
integration because of the elliptical orbit, the average bias of 
the signal in the cross gradient term is zero (see previous dis- 
cussion about Figure 15}) and the integration proceeds well. (The 
slight ridge in the back portion of the plot in Figure {17} is 
due to computer roundoff error.) 



Figure (17) - Integrated RGG Cross Gradient Output 
Three Mascons at 120 Km Spacing 
60 Km Perilune of e»0.005 Elliptical Orbit 
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The simulation was then changed so that the apolune of the 
elliptical orbit was over the three mascon set. The principal 
gradient difference output of a single RGG is plotted in Fig- 
ure {18}. Notice that the curvature of the data is in the oppo- 
site direction, and that the sensitivity and resolution of the 
data is less. This is to be expected from the increase in space- 
craft altitude over the three mascon set from 60 Km to 78 Km. 



Figure {18} - RGG Principal Gradient Difference Output 
Three Mascons at 120 Km Spacing 
78 Km Apolune of e*0.005 Elliptical Orbit 
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Figure {19} gives the integrated cross gradient output of a 
single RGG system. Again, this plot has less resolution because 
of the data integration and the higher spacecraft altitude. 



Figure {19} - Integrated RGG Cross Gradient Output 
T..ree Mascons at 120 Km Spacing 
78 Km Apolune of e=0.005 Elliptical Orbit 
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We next tried to get a feeling for the combined effects of a 
number of error sources on the performance of a single RGG sys- 
tem. We simulated an elliptical orbit with an eccentricity of 
0.005 and with the three mascon set halfway between perilune and 
apolune so that there would be a large background bias shift. We 
then introduced a periodic attitude control error perturbation 
with an amplitude of 0.1 degree and a five minute period. This 
introduced bias shift errors due to the attitude change, and ang- 
ular rate errors due to the gradient of the rate of the error. 
Figure {20} shows that the effects of the errors are significant, 
but the mascons are still resolvable. 



Figure {20} - RGG Principal Gradient Difference Output 
Three Mascons at 120 Km Spacing 
0.1 Degree Periodic Attitude Error 
e*0.005 Elliptical Orbit 
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The plot of the vertical gravity obtained from the integrat- 
ed cross gradient output of the single RGG system (See Fig- 
ure {21}) shows a considerably better picture. Although the 
resolution is still not as good because of the inherent resolu- 
tion advantage of the gradient over the gravity, there is still 
adequate resolution, and of equal importance, better signal-to- 
noise. This is partially due to the integration process which 
tends to smooth out the higher frequency components of the data 
and noise. 



Figure {21} - Integrated RGG Cross Gradient Output 
Three Mascons at 120 Km Spacing 
0.1 Degree Periodic Attitude Error 
e=0.005 Elliptical Orbit 
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However, the better signal co-noise of the integrated cross 
gradient output of the single RGG system is not just due solely 
to the integration of the data, but is due partially because the 
cross gradient output has smaller variations caused by the atti- 
tude control errors. As we can see in comparing Figure {22} with 
Figure {20}, the attitude error variations are less for the cross 
gradient output. Thus this shows that the two outputs are com- 
plementary and both should be used for best results. 



Figure {22} - RGG Cross Gradient Output 

Three Mascons at 120 Km Spacing 
0.1 Degree Periodic Attitude Error 
e*0.005 Elliptical Orbit 
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To give a dramatic presentation of the effects of very bad 
attitude control problems on the two outputs of a single RGG sys- 
tem, we increased the amplitude of the previous attitude control 
error to 1 degree. We kept the period at 5 minutes and kept the 
eccentricity of the orbit the same. Figure {23} is a plot of the 
principal gradient difference output of the single RGG system. 
One can see the mascon peaks in there, and can certainly tell 
there are three distinct peaks, but an estimate of their magni- 
tude would be of poor accuracy. 



Figure {23} - RGG Principal Gradient Difference Output 
Three Mascons at 120 Km Spacing 
1 Degree Periodic Attitude Error 
e=0.005 Elliptical Orbit 
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However, a plot of the integrated cross gradient output of 
the single RGG system (see Figure {24}), shows that this output 
is nowhere near as drastically affected, even by these very large 
error rates. 



5337-24 



Figure {24} - Integrated RGG Cross Gradient Output 
Three Mascons at 120 Km Spacing 
1 Degree Periodic Attitude Error 
e»0.005 Elliptical Orbit 
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The conclusions that can be drawn from this simulation are: 

Spacecraft altitude variations due to an elliptical orbit 
may cause some loss of sensitivity and resolution, but the bias 
shifts introduced should not prevent the extraction of the higher 
frequency gravity data. 

Spacecraft attitude tilts cause a bias shift, and should be 
rapidly recognizable and removable from the data by their appear- 
ance in the cross gradient terms. 

Periodic attitude errors caused by malfunction of the atti- 
tude control system will cause significant data reduction errors 
if the attitude error amplitudes are large (greater than 1 
degree) and of rapid period (less than 5 minutes). (Since they 
are time varying, they will show up as angular tilts in the cross 
gradient outputs and as angular rate gradients of the same peri- 
odicity but quadrature phase in the trace.) 

A single RGG sensor oriented along the orbital axis to mea- 
sure the gravity gradients in the plane of the orbit can use the 
lunar orbit track to track spacing of about 30 Km to obtain good 
resolution of mascons in both the along track and cross track 
directions. The cross gradient output of the RGG should be inte- 
grated with the spacecraft velocity to obtain a contour map of 
the vertical gravity field, and then the higher resolution (but 
slightly distorted) principal gradient difference output of the 
RGG should be used to enhance the resolution of the data set for 
geophysical interpretation purposes. 
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